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1  his  study  i mitigates  the  polarimetrir  behaviour  of  canonicals  w;thm  •!:!;> 
m-o’mtion  ceils.  1  he  aim  of  this  investigation  is  to  ascertain  whet  Imr  the  c. 

•  c  ,1  resolution  cell  can  be  determined  by  examination  oi  the  polarma ■'  i  ,c  •  ; 

canonicals  addressed  in  this  t hesis  effort  are  the  dihedral  and  the  iiatpiato.  , : . ’ ; . ■ 

'  associated  theory  and  software  is  readily  applicable  to  other  j ■  «. 

The  main  software  package  used  during  this  thesis  is  Sarloor".  alt  lacath 
• : : ’  tonal  ancillary  software  is  presented  such  that  the  output  of  Sar  Tool'"*  ■  m. 

: e; u '‘sriit ei  1  cm  the  Poincare  sphere. 

Finally,  both  single  level  and  multi-level  percept ron  networks  are  used  to  pm 
<  e>.<  the  resolution  cell  output  and  return  a  decision  on  the  contents  ol  the  ceil.  1 
preliminary  study  of  the  polarimetric  and  percept  ron  combination  provides  so::.'' 
valuable  insight  into  the  possible  potential  of  polarirnetrics  and  Artificial  Intelligence 
in  the  area  of  target  identification. 
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TARGET  CLASSIFICATION  USING  SYNTHETIC  APERTURE 

RADAR 

POLARIMETRIO  DATA 


I.  Introduction 


I.l  Background 

The  origins  of  radar  theory  ran  be  traced  to  the  early  part  of  the  cc'ntury: 
however,  practical  radar  systems  only  gained  prominence  during  World  War  II  where 
they  were  used  extensively  by  the  British  Chief  of  Naval  Operations  as  Earlv  Warning 
(LAV.)  devices.  Even  in  those  early  days  operators  recognized  that  the  theory  of 
radar  was  not  an  exact  science  [ij.  Variations  in  target  detection  ranges  and  in  false 
alarm  rates  led  to  a  flurry  of  postwar  studies  which  ultimately  resulted  in  the  first 
generation  of  radar  systems.  These  radars  are  commonly  referred  to  as  conventional 
radars. 

1.1.1  Conventional  Radars.  Chapter  2  of  this  document  refers  the  interested 
reader  to  texts  which  review  the  theory  of  conventional  radars.  From  a  systems 
point  of  view  the  limitation  of  conventional  radars  is  their  inability  to  identify  cer¬ 
tain  tactical  targets.  This  lack  has  limited  the  use  of  conventional  radars  to  that  of 
target  detection  systems.  Simple  target  detection  has  proven  less  than  satisfactory 
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to  both  military  and  civilian  users  of  radars.  From  ?  military  perspective  modern 
weapons  have  resulted  in  increased  Missile  Engagement  Zones  (MF.Z)  and  the  sub¬ 
sequent  development  of  "stand  off  and  attack"  doctrines  which  necessitate  brvond 
visual  range  identification  of  potential  targets.  Similarly,  the  increasing  use  cf  space 
by  civilian  agencies,  such  as  Land  Satellite  (LANDSAT)  has  necessitated  the  long 
range  identification  of  certain  Geophysical  features,  for  example,  coastlines,  weather 
patterns,  and  crops. 

1.1.2  SAR  Systems.  Since  the  mid  19G0s  the  limitat;ons  of  conventional 
radars  and  the  advent,  of  new  technologies  has  prompted  considerable  interest  in  the 
development  of  SAR  systems.  Several  good  texts  which  review  the  theory  of  SAR 
systems  are  referenced  in  the  literature  review  section  of  this  paper.  Bas'cally  SAR 
systems  achieve  better  resolution  performance  than  conventional  radars  due  ‘o  the 
formers  ability  to  use  the  platform's  path  as  if  it  were  a  physical  aperture.  The  use 
of  the  platform’s  flight  path  as  a  synthetic  aperture,  coupled  with  advanced  process¬ 
ing  algorithms  and  sizeable  computer  resources,  results  in  resolutions  significantly 
better  than  those  attainable  by  conventional  radars.  Consequently  SAR  systems 
have  become  the  “second  generation  of  radars”  and  have  been  widely  accepted  for 
use  in  both  military  and  civilian  applications. 

1.1.3  Polarimetric  Research.  Despite  their  improved  performance  over  con¬ 
ventional  radars,  modern  SAR  systems  are  not  able  to  perform  certain  tactical  tasks, 
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for  example,  the  identification  of  tanks  and  aircraft.  Current  research  indicates  that 
the  only  way  SAR  systems  will  he  able  to  achieve  the  performance  necessary  for  these 
tasks  is  by  exploiting  the  polarization  information  contained  in  the  radar  return. 

During  the  early  days  of  SAR  development,  scientists  noted  that  the  returns 
from  targets  were  invariably  de-polarized  (i.e.  they  returned  to  the  receiver  at  a 
ditferent  polarization  to  that  which  they  were  transmitted  on).  Further  investigation 
revealed  that  this  process  of  de-polarization  was  not  a  completly  random  mechanism. 

This  revelation  raised  questions  as  to  whether  or  not  greater  target  identifica¬ 
tion  (and  hence  classification)  could  be  attained  if  the  polarization  of  radar  returns 
was  included  in  the  process  of  target  analysis.  Consequently,  the  study  of  Polarimet- 
rics  began.  According  to  radar  historians,  this  research  started  in  the  early  1950's 
and  continued  into  the  1960’s  when  it  lapsed  due  to  the  following  reasons: 

l)Polarization-target  dependence  was  an  incompletely  understood  phe¬ 
nomenon;  2)  Because  polarization  diversity  techniques  require  antenna 
polarization  control  or  agility  on  transmit  and/or  dual  polarization  chan¬ 
nels  on  reception,  the  radar  system  complexity  increased  significantly  and 
3)  the  knowledge  of  polarization  behaviour  of  complex  target  requires  ex¬ 
tensive  and  expensive  dual  polarization  measurements  [2:42]. 

In  1970,  Iluynen  [3]  revived  interest  in  polarimetric  studies.  This  revived 
interest,  coupled  with  important  theoretical  and  technical  advances  in  polarimetric 
analysis,  continues  today.  The  ultimate  objective  of  this  research  is  to  significantly 
improve  the  classification  capabilities  of  SAR  systems. 
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1.1.4  Problem  Statement.  The  objective  of  this  study  is  to  analyze  the  po- 
larimetric  behaviour  of  certain  canonical  shapes  within  a  single  resolution  cell  of  a 
SAR. 


1.1.5  Scope.  The  study  of  polarimetrics  is  extensive  and  encompasses  several 
mathematical  and  scientific  disciplines  pertinent  to  target  definition.  However,  due 
to  the  time  constraints  imposed  by  the  thesis  program,  the  scope  of  this  thesis  effort 
will  be  limited  to: 

1.  Establishing  the  accuracy  of  the  SarTool£m  software  in  predicting  the  polari- 
metric  behaviour  of  the  canonical  shapes  used  within  this  thesis. 

2.  Investigating  ways  in  which  SarTool£m  data  can  be  effectively  displayed  using 
either  the  Poincare  sphere  or  the  Matrix  Laboratory  Software  (MATLAB(m) 
or  a  combination  of  the  two. 

3.  Examining  the  polarimetric  behaviour  of  individual  canonical  shapes  within  a 
single  SAR  resolution  cell. 

4.  Examining  the  behaviour  of  two  or  more  canonical  shapes  co-located  within  a 
single  SAR  resolution  cell. 

5.  Investigating  mean  for  discriminating  between  various  shapes  and  combinations 
of  shapes  within  a  resolution  cell. 
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1.2  Summary 

Since  the  development  of  the  first  radar  system  in  World  War  II  there  has  been 
much  ongoing  research  directed  at  improving  the  resolution  performance  of  opera¬ 
tional  radar  systems.  This  research  is  of  paramount  importance  as  the  resolution 
performance  of  a  radar  to  a  large  extent  determines  the  overall  capability  of  that 
radar. 

Current  SAR  systems  have  resolution  capabilities  far  in  excess  of  conventional 
radars.  Nevertheless,  even  with  this  enhanced  resolution,  today's  SAR  systems  are 
incapable  of  performing  several  critical  functions,  for  example,  the  identification 
of  tactical  targets.  Furthermore,  there  are  physical  and  system  constraints  which 
preclude  these  radars  from  attaining  greater  target  classification  performance  unless 
polarimetrics  are  incorporated  into  the  target  discrimination  process. 

Unfortunately,  the  polarimetric  behaviour  of  targets  is  an  incompletely  under¬ 
stood  phenomenon.  However,  there  is  sufficient  understanding  of  polarimetrics  to 
indicate  that  this  avenue  of  research  may  ultimately  lead  to  target  classification  per¬ 
formances  far  in  excess  of  those  currently  attainable.  Such  an  improvement  in  radar 
performance  would  significantly  increase  the  usefulness  of  SAR  systems  to  civilian 
and  military  users  alike. 

This  thesis  effort  is  directed  at  increasing  the  knowledge  of  polarimetrics  by 
investigating  the  polarimetric  behaviour  of  certain  canonical  shapes  within  a  sin¬ 
gle  SAR  resolution  cell  and  by  examining  ways  in  which  the  resulting  polarimetric 
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information  can  be  used  in  the  decision  process. 


II.  Literature  Search 


1  Introduction 

As  mentioned  in  Chapter  1,  the  study  of  polarimetrics  encompasses  a  large 
number  of  research  activities;  accordingly,  before  progressing  further,  a  literature 
review  is  necessary  to  identify  the  material  relevant  to  this  thesis  effort. 

2. 2  Scope 

The  literature  reviewed  here  provides  an  overview  of  the  development  and  po¬ 
tential  of  Polarimetric  studies  as  they  pertain  to  the  current  areas  of  Radar  Imaging. 

2.3  Organization 

To  fully  appreciate  the  relevance  and  theory  of  Polarimetrics,  one  must  first 
understand  how  the  study  of  Polarimetrics  applies  to  Radar  imaging.  Moreover, 
as  the  study  of  Polarimetrics  is  primarily  theoretical  (with  only  a  few  practical 
examples),  the  preponderance  of  papers  on  the  subject  deal  mainly  with  the  math¬ 
ematical  modelling  of  expected  responses.  Consequently,  there  is  a  need  to  review 
the  mathematical  tools  and  modelling  processes  currently  being  proposed  by  the 
leading  authorities  in  this  area  of  study.  Accordingly,  this  Literature  Review  will  be 
organized  in  the  following  manner: 

1 .  Background.  This  section  of  the  review  will  consist  of  the  following  subsections: 
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(a)  Conventional  Radars. 


(b)  SAR  Systems. 

2.  Mathematical  Tools.  This  section  will  briefly  introduce  the  literature  regarding 
the  currently  accepted  mathematical  tools  of  Polarimetrics,  in  particular  the 
following  topics  will  be  reviewed: 

(a)  The  Polarization  Ellipse. 

(b)  The  Poincare  Sphere. 

(c)  The  Scattering  Matrix. 

3.  Results  of  Previous  Polarimetric  Studies.  This  section  will  provide  a  brief 
overview  of  several  technical  papers  which  compare  actual  Radar  Images  with 
those  obtained  via  modelling  with  the  Mathematical  tools  mentioned  above. 

2.4  Background 

2.4-1  Conventional  Radars.  The  theory  of  conventional  radars  is  well  known 
and  is  the  subject  of  several  textbooks.  A  particularly  excellent  treatment  of  the 
theory  of  conventional  radars  is  found  in  [4]. 

As  mentioned  in  Chapter  1,  the  salient  limitation  of  conventional  radars  is 
their  lack  of  both  down  range  and  cross  range  resolution.  For  example,  typical 
radar  parameters  for  an  aircraft  search  radar  are:  1)  RF  3  GHZ,  2)Pulse  Width  .5 
microseconds,  3)  Aperture  size  5  meters  and  4)  “look”  angle  of  3  degrees.  Application 
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of  conventional  radar  theory  via  equations  (2.1)  and  (2.2)  show  that  the  best  possible 


down  range  and  cross  range  resolutions  (on  a  target  at  100km)  are  75  meters  arid 
610  meters  respectively  [5]: 


A  Rr 


TC 

2  cos  a 


•i  i  i 


where, 


Rr  =  down  range  resolution 
t  =  pulse  width 
a  =  look  angle 


ARc  =  RjS4  (2.2) 

where, 

Rc  =  cross  range  resolution 
A  =  radar  wavelength 
L  =  length  of  flight  path 

2.4-2  SAR  Systems.  The  above  resolutions  are  capable  of  providing  simple 
functions,  for  example,  target  detection  or  coastline  mapping.  However,  more  ad¬ 
vanced  capabilities,  such  as  target  and  crop  identification,  require  considerably  better 
resolutions  than  those  attainable  via  conventional  radars.  Accordingly  OAK  systems 


were  developed  which  which  could  yield  resolutions  less  than  1  meter  [6].  Chapter 
1  points  out  that  these  much  improved  resolutions  are  possible  due  to  significant 
technological  advances  and  improvements  in  processing  techniques.  A  definitive 
treatment  of  these  techniques  is  beyond  the  scope  of  this  thesis:  however,  a  brief 
review  of  SAR  theory  is  included  in  the  following  chapter.  A  theoretical  aspect  of 
immediate  relevance  is  that  modern  SAR  systems  exploit  all  available  radar  param¬ 
eters  except  one.  The  exploited  parameters  are:  time,  frequency,  phase,  hewing  and 
amplitude.  The  one  parameter  not  used  by  current  SAR  systems  is  the  polarization 
of  the  Electro-Magnetic  (EM)  wave. 

ine  need  for  even  better  resolution  performance  than  is  achievable  with  con¬ 
ventional  SAR  systems  has  led  to  the  study  of  polarimetrics  which,  as  mentioned 
earlier,  involves  understanding  the  polarization  behaviour  of  E.M.  radiation.  An 
integral  part  of  this  essentially  theoretical  study  is  the  use  of  several  mathematical 
tools. 


2.5  Mathematical  Tools. 

Much  of  the  research  into  polarimetrics  is  theoretical  and  based  on  the  de¬ 
velopment  and  use  of  mathematical  tools.  Accordingly,  this  section  will  review  the 
literature  concerning  these  tools. 

Many  authors  such  as  [7]  and  [8],  have  spent  considerable  research  effort  in 
using  and  improving  the  math  and  processing  techniques  associated  with  polarimet- 
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tic  analy  The  consensus  among  these  authors  is  that  the  following  three  basic 
mathematical  tools  are  fundamental  to  the  study  of  polarimetrics: 


1.  The  Polarization  Ellipse. 

2.  The  Poincare  Sphere 

3.  The  Scattering  Matrix 

2.5.1  Polarization  Ellipse.  All  polarization  states  can  be  described  using  the 
Polarization  Ellipse  [9],  an  example  of  which  appears  in  Figure  2.1.  A  brief  descrip¬ 
tion  of  the  significance  of  Figure  2.1  is  given  below: 

1.  Ellipticity.  The  ellipticity  angle  \  reveals  the  polarization  angle  of  the  wave. 
For  example,  an  ellipticity  angle  of  0°  would  denote  linear  polarization  whereas 
an  angle  of  45°  would  indicate  circular  polarization.  Moreover,  the  sign  of  the 
ellipticity  angle  denotes  the  “handedness”  of  the  wave  i.e.  whether  the  wave 
is  travelling  towards  an  observer  (right  handed)  or  away  from  an  observer  (left 
handed). 

2.  Ellipse  Orientation  Angle.  The  Ellipse  Orientation  Angle  'I'  represents  the 
“tilt.”  or  orientation  of  the  ellipse  relative  to  a  defined  axis.  The  tilt  or  orien¬ 
tation  of  an  EM  wave  is  an  important  consideration  in  defining  a  wave  in  a 
specified  coordinate  system. 


2-5 


h 


Figure  2.1.  The  Polarization  Ellipse  [2 


2.5.2  Poincare  Sphere.  The  Poincare  Sphere  is  described  in  considerable  de¬ 
tail  in  most  papers  pertaining  to  polarimetric  studies.  Kennaugh  is  generally  ac¬ 
cepted  as  being  the  leading  authority  on  the  Poincare  sphere,  and  during  last  30 
years  he  has  written  some  excellent  papers  on  its  applications  to  polarimetrics  [10], 
The  mapping  of  the  Polarization  ellipse  onto  the  Poincare  sphere  is  accomplished 
as  shown  in  Figure  2.2.  Basically,  each  ellipticity  angle  \  and  orientation  angle 
$  can  be  uniquely  mapped  to  a  point  on  the  sphere  having  longitude  2$  and 
latitude  2\  .  Examination  of  the  sphere  reveals  that  the  poles  represent  the  circular 
polarizations,  the  equator  represents  the  linear  polarizations  and  the  hemispheres 
represent  the  “handedness”  of  the  wave.  The  power  in  the  EM  wave  is  proportional 
to  the  radius  of  the  sphere.  One  of  the  great  strengths  of  the  Poincare  sphere  is  that 
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Figure  2.2.  The  Poincare  Sphere  [8] 


it  can  be  used  to  show  the  statistical  or  stochastic  variations  of  polarizations.  For 
example  the  shaded  area  in  Figure  2.2  indicates  the  degree  of  polarization  jitter  in 
the  measured  wave’s  ellipticity  and  orientation  angles. 

2.5.3  Scattering  Matrix.  There  is  a  considerable  amount  of  scientific  liter¬ 
ature  devoted  to  the  scattering  matrix  because  mathematical  techniques,  coupled 
with  new  processing  systems,  allows  the  scattering  behaviour  of  any  object  to  be 
modelled  as  a  complex  matrix.  The  theory  is  that  if  the  amplitude  and  polarization 
of  an  EM  wave  is  measured  (before  and  after  it  is  scattered  by  a  target)  then  a 
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matrix  can  be  formulated  which  can  simulate  the  same  behaviour  as  the  scat  ten  r. 
The  logical  extrapolation  of  this  is  that  if  a  "data  bank"  of  these  matrices  exists, 
each  identified  with  particular  canonical  shape,  then  any  composite  shape  can  be 
identified  by  a  process  of  matrix  comparison.  Equally  important,  if  the  scattering 
behaviour  of  basic  shapes  is  known  (particularly  in  a  polarization  sense),  then  ' u,, 
adjacent  objects  may  be  distinguished  from  each  other  by  making  finite  comparisons 
ut  their  scattering  behaviour  !ll]- 

A  particularly  useful  ilhist  ration  showing  the  applicability  of  matrix  mat 'ne¬ 
matics  to  scattering  problems  can  be  found  in  a  1983  article  by  Becker  IT.  Fur¬ 
thermore,  a  paper  written  by  Polikarpov  [13]  in  1985  shows  how  the  seatteiing 
matrix  can  be  measured  using  linear  regression  techniques  and  the  method  of  least 
-quares.  The  following  chapter  will  address  portions  of  this  theory  and  will  >!mw 
how  matrix  mathematics  can  he  combined  with  the  Poincare  sphere  to  assist  in  the 
analysis  of  polarimetric  data. 

2.6  Results  of  Previous  Polarimetric  Studies 

Several  papers  have  analyzed  the  interrelationships  between  a  radar's  polar¬ 
ization  and  the  associated  return  from  a  target.  The  following  section  summarizes 
the  three  most  topical  of  those  papers. 

2.6.1  Radar  Backscattering  From  Trees.  This  paper  mentioned  earlier  that 
both  civilian  and  military  authorities  were  involved  in  SAR  systems  and  research 
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into  polarimetrics.  A  paper  which  typifies  the  interest  of  the  civilian  -cctor  w.i- 
written  in  1988  as  a  result  of  Polarimetric  studies  conducted  on  various  tree  tcpe- 
(hiring  the  1987  growing  season  [14].  This  study  concluded  that  the  measured  Radar 
Cross  Section  (RCS)  of  trees  varied  according  to  the  selected  radar  polari/at  ion. 
Moreover,  by  a  series  of  exhaustive  measurements  the  author  was  able  to  diow  t  hat 
discernible  difference  existed  between  the  RCS  behaviour  of  different  species  of  tic- 
These  conclusions  are  significant  as  they  demonstrate  that  the  use  of  polarimet  rn 
data  may  lead  to  classification  capabilities  well  beyond  those  attainable  bv  even  ’be 
highest  resolution  non  polarimetric  SAR  systems. 

2.6.2  Dihedral  Corner  Reflectors.  A  seemingly  inordinate  amount  of  the  un¬ 
classified  polarimetric  data  is  concerned  with  studies  which  examine  the  RCS  chnrae- 
teristics  of  dihedral  corner  reflectors  [15].  However,  the  motive  behind  1 1 . i -  apparent 
pre-occupation  with  dihedraia  is  that  they  typify  man  made  structures.  Such  •'true- 
lures  tend  to  be  angular  and  have  high  RCS's  whereas  natural  objects  usually  have 
"softer  ”  lines  and  lower  RCS's.  This  important  Geometrical  observation  is  instru 
mental  in  the  matrix  formulations  mentioned  earlier  in  this  paper.  Moreover,  at  a 
gross  level  this  differing  RCS  behaviour  explains  why  many  military  platforms  e.u. 
Stealth  Aircraft  have  predominantly  curved  shapes  [16], 

Two  very  good  discussions  on  the  importance  of  dihedrals  to  polarimetric  stud¬ 
ies  are  to  be  found  in  [17]  and  [18].  fn  the  former  paper  the  author  demonstmtes 
excellent  correlation  between  the  predicted  and  the  measured  reaction  of  a  dihedral 


corner  reflector  to  varying  degrees  of  polarization.  This  encouraging  result  is  some¬ 
what  tempered  by  [19]  which  concludes  that  although  math  models  can  predict 
polari metric  behaviour,  those  same  models  can  give  misleading  information  if  the 
target  geometry  is  not  carefully  considered.  This  “constrained  ’uccess"  typifies  the 
current  status  of  Polarimetric  studies. 

3.7  Summary 

The  study  of  Polarimetric  radar  theory  is  an  emerging  field  with  its  own  subset 
of  scientific  jargon,  an  understanding  of  which  is  necessary  in  order  to  comprehend 
the  books  and  journals  written  on  this  subject.  Similarly,  an  appreciation  of  the  dif¬ 
ferences  in  performance  between  conventional  radars  and  SAR  systems  is  central  to 
understanding  the  transition  between  the  well  known  technology  of  the  conventional 
systems  and  the  “still  being  investigated”  technology  of  the  Polarimetric  SAR  sys¬ 
tems.  Finally,  there  are  several  mathematical  tools,  in  varying  stages  of  development  , 
which  can  be  used  to  predict  the  scattering  behaviour  of  canonical  shapes. 

Current  attempts  to  synthesize  the  tools  and  the  technology  into  a  reliable 
predictive  mechanism  have  met  with  guarded  success.  Complete  success  will  be 
difficult  to  achieve  due  to  the  many  parameters  involved  and  the  statistical  nature 
of  EM  propagation.  Nevertheless,  the  latent  benefits  of  integrating  Polarimetrics 
into  the  radar  recognition  and  identification  process  are  considerable  and  have  the 
potential  to  revolutionise  the  performance  capabilities  of  future  SAR  systems. 


2-10 


III.  Theory 


3.1  Introduction 

As  mentioned  in  Chapter  2  there  is  a  need  to  present  the  applicable  SAR  theory 
so  that  the  methodology  and  results  of  this  thesis  effort  can  be  placed  in  the  context  of 
their  potential  to  contribute  to  the  overall  improvement  of  SAR  performance.  There 
is  also  a  need  to  address  how  practical  systems  could  use  the  derived  polarimetric 
data  in  the  decision  or  classification  process.  Accordingly, this  chapter  will  consist  of 
the  following  sections: 

1.  Non  Polarimetric  SAR  Resolution  Limits. 

2.  Thesis  Methodology. 

3.  Perceptron  Network  Theory. 

3.2  Non  Polarimetric  SAR  Resolution  Limits 

3.2.1  Categories  of  SAR  Systems.  SAR  systems  can  be  classified  into  one  of 
three  major  categories.  These  categories,  and  a  brief  discussion  of  their  operational 
concepts,  are  as  follows: 

1.  Strip  Map  SAR.  In  strip  map  SAR  the  radar’s  real  beam  is  fixed  in  both 
azimuth  and  elevation.  Accordingly,  the  area  to  be  mapped  is  a  function  of 
the  platform’s  flight  path. 
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2.  Spotlight  SAR.  In  spotlight  SAR  the  radar  is  steered  such  that  the  area  of 
interest  is  illuminated  throughout  the  processing  interval. 

3.  Scan  Mode.  Scan  mode  SAR  is  a  mixture  of  both  strip  and  spotlight  SAR.  In 
scan  mode  the  area  to  be  illuminated  is  a  function  of  both  platform  path  and 
beam  steering. 

4.  Inverse  SAR  (ISAR).  The  theory  associated  with  ISAR  has  much  in  common 
with  the  previously  mentioned  SAR  modes.  The  single  most  fundamental  dif¬ 
ference  between  SAR  and  ISAR  is  that  the  latter  relies  on  the  target’s  motion, 
rather  than  the  platform’s  motion,  to  perform  the  mapping  function. 

The  radar  prediction  software  used  for  this  thesis  effort  is  the  Synthetic  aper¬ 
ture  radar  prediction  Tool  (SarToolfm).  SarTooltm  was  developed  by  The  Analytic 
Sciences  Corporation  (TASC)  and  is  capable  of  simulating  all  of  the  above  radar 
modes.  A  brief  description  of  the  salient  details  of  SarTooltm  is  contained  in  Ap¬ 
pendix  A.  For  the  purposes  of  this  thesis,  the  software  input  files  were  configured 
to  simulate  the  ISAR  mode  of  operation.  This  involved  holding  the  SAR  platform 
stationary  and  spotlighting  a  rotating  pedestle  upon  which  the  targets  were  located. 
Accordingly,  the  following  theory  will  address  the  resolution  cell  sizes  applicable  to 
spotlight/ISAR  systems  and  the  sampling  requirements  of  those  cells. 

3.2.2  Cross  Range  Resolution.  The  cross  range  resolution  of  a  SAR  system 
can  be  explained  in  terms  of  physical  constraints  (i.e.  antenna  size)  or  in  terms 
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Point  Targets 


v  =  —uy 

Figure  3.1.  Doppler  Cross  Range  Resolution [20] 

of  discernible  differences  in  target  doppler.  The  latter  approach  is  more  applicable 
to  ISAR:  accordingly,  the  following  development,  [20]  defines  cross  range  resolution 
primarily  in  terms  of  target  doppler.  Wehner’s  doppler  development  centers  on 
Figure  3.1,  where: 

R  =  range  to  centroid  of  targets 
Vp  =  straight  line  velocity  of  platform 
u>  =  angular  velocity  of  platform 
'F  =  equivalent  rectangular  beamwidth 

In  accordance  with  Figure  3.1  target  1  and  target  2  will  have  respective  doppler 
shifts  of  —(2 f/c)wy  and  -f-(2 f/c)wy,  where  /  is  the  radar’s  center  frequency.  Ac- 


cordinglv,  the  difference  in  doppler  frequency  between  the  two  targets  is  given  by 


equation  (3.1)  shown  below: 


°fd  =  -uyf  (3.1) 

c 

Rearranging  equation  (3.1)  and  substituting  A f  =  afd  —  system  doppler  fre¬ 
quency  resolution  leads  to  the  cross  range  resolution  (Arc)  being  defined  by  the 
following  equation: 


Arc 


c 

S Td 


(3.2) 


Recognizing  that  for  ideal  processing  of  a  uniform  signal  6fd  =  l/T  then  equa¬ 
tion  (3.2)  can  be  re-cast  in  any  one  of  the  following  forms: 


Arc 
A  rc 


Arc 


A 

2 ujT 

A 

2A 

A 

2A$ 


(3.3) 

(3.4) 


(3.5) 


where 

A3>  =  look  angle  =  u>T 

Wehner  also  shows,  via  integration  of  the  response  for  a  rectangular  beam 


antenna,  that  the  cross  range  resolution  can  also  be  expressed  by  equation  (3.6) 
which  is  shown  below: 


(3.6) 

where 

l  =  aperture  length  of  the  physical  antenna  in  the  azimuth  plane 

3.2.3  < lant  Range  Resolution.  SAR  systems  rely  on  some  form  of  transmitter 

pulse  coding  to  achieve  slant  range  resolution  improvement  over  conventional  radars. 
In  SarTool(m  this  pulse  coding  takes  the  form  of  CHIRP  or  pulse  compression  coding. 
The  frequency  -  time  representation  of  a  chirp  pulse  is  shown  in  Figure  3.2 

As  illustrated  in  Figure  3.2  the  chirp  pulse  has  a  unique  characteristic  over 
every  finite  time  interval.  If  this  were  not  the  case  then  SAR  radars  could  only  achieve 
slant  range  resolutions  of  CTi/2  where  T\  =  transmitter  pulse  width.  However,  to 
satisfy  the  Nyquist  sampling  theorem  each  echo  from  the  linear  FM  pulse  must  be 
sampled  at  /3  times  per  second.  This  leads  to  the  following  simple  expression  for 
slant  range  resolution: 


(3.7) 
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Figure  3.2.  Chirp  or  Pulse  Compression 


where 

=  slant  range  resolution 

The  actual  sampling  criteria  for  both  cross  and  slant  range  resolution  can  be 
derived  using  Figure  3.3.  To  adequately  sample  in  slant  and  range  the  following 
relationships  must  apply: 


A  ^  CT 1 

'.Ar,>  2 

(3.8) 

?7cArc  >  Rrp 

(3.9) 

where 
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i)f  =  samples  collected  during  T\ 

ijc  =  samples  collected  along  r  ip 
Rearranging  equation  (  3.8)  leads  to 


7*  > 


CTX 

2  Ar, 


=  7\A 


where 


13.10) 


A  =  chirp  pulse  bandwidth 
Similarly,  re-arranging  equation(3.9)leads  to: 


7c  > 


Rip 

Ary 


(3-11) 


or 


7c  > 


2R&(prpe 

A 


(3.12) 


The  above  discussion  enables  the  maximum  resolution  capabilities  of  a  conven¬ 
tional  SAR  system  to  be  calculated  given  that  the  necessary  sampling  criteria  are 
adhered  to.  Both  equation  (3.10)  and  equation(3.12)  are  predicated  on  one  complex 
sample  per  resolution  cell.  SarTooltm  satisfies  the  number  of  transmitted  pulses  and 
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Figure  3.3.  SAR  Data  Sampling  [20] 
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the  necessary  increment  angles  by  equation  (3.13)  and  equation  (3.14). 


np 
A  0 


ARC 

A0 

C 

2  nadacos{$) 


(3.13) 

(3.14) 


where 


np  =  number  of  transmitted  pulses 
ARC  =  desired  portion  of  angular  arc  of  the  turntable 
A0  =  angular  increment  in  aspect  angle 

na  =  number  of  azimuth  bins 

da  =  azimuth  resolution 

$  =  elevation  boresight 

Of  note  is  that  although  the  preceeding  discussion  determines  the  number  of 
pulses  required  for  a  particular  resolution  it  does  not  address  the  use  of  polarimetric 
information.  Accordingly,  the  resulting  resolution  relationships  are  absolute  limita¬ 
tions  of  conventional  SAR  systems.  As  mentioned  in  Chapter  1  this  thesis  effort 
is  directed  at  ascertaining  whether  or  not  intra  resolution  cell  classification  can  be 
achieved  via  the  use  of  polarimetric  data.  Consequently  the  following  section  will 
address  the  methodology  used  to  collect  and  analyze  the  intra  resolution  cell  polari¬ 
metric  data. 
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3.3  Thesis  Methodology 

3.3.1  Introduction.  As  indicated  in  the  previous  section  of  this  chapter,  cur¬ 
rent  SAR  systems  have  target  identification  limits  that  are  unassailable  without  the 
use  of  polarimetrics.  Consequently,  this  thesis,  (which  is  intended  to  he  the  first 
of  several  similar  efforts)  is  directed  towards  assessing  the  viability  of  including  po¬ 
larisation  within  classification  algorithms.  In  particular  this  thesis  addresses  the 
following  areas  of  concern: 

1.  The  validation  of  SarTool<m. 

2.  The  development  of  a  composite  software  package  to  aid  subsequent  investiga¬ 
tions  into  polarimetric  behaviour. 

3.  The  examination  of  polarimetric  data  from  single  and  double  bounce  canonicals 
within  a  single  resolution  cell. 

3.3.2  The  Validation  of  SarTool  tm.  The  sponsor  of  this  thesis  was  inter¬ 
ested  in  determining  how  accurately  the  documentation  of  SarTooltm  reflected  the 
performance  of  the  software.  Due  to  the  limited  time  available,  and  the  secondary 
nature  of  this  task,  it  was  agreed  to  validate  only  the  dihedral  section  of  the  doc¬ 
umentation  and  software.  This  validation  entailed  programing  MATHCAD(m  and 
Mathematical  (both  of  which  are  math  packages)  with  the  equations  purportedly 
embedded  in  the  SarTooltm  algoiithms.  To  facilitate  this  validation,  the  parameters 
used  in  the  math  packages  and  SarTool*m  corresponded  to  those  use  in  an  IEEE 


article  written  by  Griessner  and  Balanis  [17]  and  the  azimuthal  scope  was  limited  to 
90  degrees.  The  actual  comparison  was  achieved  by  •‘shipping"  output  SarTool,m 
data  and  plotting  that  data  within  a  Mathematicaim  routine. 

In  addition  to  the  above  validation  several  experiments  were  conducted  in 
order  to  ascertain  how  responsive  SarTool'"1  was  to  varying  input  configurations. 
The  resulting  plots  and  associated  findings  of  the  above  investigations  are  contained 
in  Chapter  IV  of  this  report. 

3.3.3  .4  Composite  SarTool  tm  Software  Package.  In  the  early  stages  of  this 
thesis  it  became  evident  that  SarTool(m,  in  its  pristine  form,  would  not  readily 
support  an  investigation  . polarimetrics.  The  major  shortfalls  of  SarTool!m  were 
assessed  as  involving  the  input  reflector  files  and  the  output  image  files.  The  input 
reflector  files  are  the  means  whereby  the  canonicals  are  defined  and  then  moved 
and  rotated  within  the  radar’s  field  of  view.  These  files  are  normally  formed  as 
outputs  of  the  EUCLID  solid  modeller  sorfware  which  was  unavailable  to  the  sponsor 
and  attempts  to  alter  the  reflector  files  dynamically  v.itlicut  the  Euclidean  software 
proved  to  be  impractical.  The  normal  output  of  SarToolim  consists  of  parametric 
tables  and  image  files,  neither  of  which  were  in  a  suitable  form  for  polarimetric 
analysis.  As  a  result  of  these  shortcomings  the  software  packages  described  below 
were  designed  and  built. 


3-11 


1.  Reflector  File  Software.  This  software  enables  the  operator  to  select  any  com¬ 
bination  of  SarTooF"1  canonical  shapes  and  translate  or  rotate  those  shapes 
within  the  resolution  cell.  The  reflector  file  software  can  then  be  used  as  an 
input  to  a  commercial  software  package  called  “did”  which  allows  visual  in¬ 
spection  of  the  canonicals. 

2.  Shell  program.  The  shell  program  is  a  routine  which  converts  the  above  ASCI 
reflector  files  into  binary  files  and  then  'Tuns"  SarTooF171  on  each  of  those 
files.  The  resulting  output  SarToolim  data  files  are  then  "stripped”  and  moved 
by  File  Transfer  Protocol  (FTP)  from  the  SUN  workstation  to  a  Macintosh 
computer  where  they  are  processed  by  the  Mathematical  routine. 

3.  Strip  Software.  The  software  which  performs  the  "stripping"  function  is  run 
from  within  the  shell  programme  and  is  used  to  transform  the  output  from 
the  SarTooF771  Range  Profile  Data  Base  (RPDB)  files  into  the  'fin  phase"  and 
"quadrature”  components.  These  in  phase  and  quadrature  files  are  subse¬ 
quently  used  as  the  input  to  the  Mathem?ticatm  routine. 

4.  Mathematica*771.  A  routine  was  written  within  Mathematica'771  which  trans¬ 
formed  the  stripped  SarTooFm  data  into  polarization  plots  and  as  "animated” 
or  stationary  polarization  points  on  the  Poincare  sphere.  The  Maihematica*771 
routine  was  based  on  Boerner’s  article  [21]  and  is  further  described  in  the 
following  paragraphs. 
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3.3.4  The  Poincare  Sphere.  The  Poincare  sphere  was  hrietiv  int roduced  iu 
Chapter  2  of  this  report  and  was  identified  as  being  of  considerable  use  in  t  lie  studv 
of  polarimetrics.  The  applicability  of  the  Poincare  sphere  to  the  study  of  polarization 
is  due  to  its  ability  to  completely  characterize  the  polarimetric  behaviour  of  a  radar 
target.  This  is  achieved  by  plotting  the  optimal  polarizations  (co- polarization  (CO- 
POL)  and  cross  polarization  (X-POL))  on  the  Poincare  sphere.  Kennaugh  was  :  he 
first  to  identify  these  "optimal  polarizations"  [22]  and  his  results  are  used  in  the 
development  of  Boerner’s  article  the  relevant  theory  of  which  is  shown  below: 

Assuming  that  no  losses  are  incurred  then  the  transformation  from  one  orthog¬ 
onal  polarization  to  another  may  be  expressed  as  shown  by  equation! 3. Is)  below: 


[S'(  A'B')]  =  [r]T[S(.4fl)][I'J 


o.r.i 


where: 


T  is  the  normalized  unitary  transformation  matrix  which  can  be  ex¬ 
pressed  in  the  form  of  equation(3. 16): 


[T]  =  (i  +  pp*)-|/2 


1  -p* 

p  1 


(3.16) 


p  =  the  complex  transformation  parameter. 
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Substituting  equation(3. 16)  into  equation(3. 15)  and  expanding  leads  to  the 


following  equations(3. 17,  3.18,  3.19,  3.20) 


S-' 

JA‘A‘ 

=  (1+PPT 

1  [-5.4.4  +  P2Sbb  +  P(  SaB  +  S/3  4 )] 

(3.17 

S' 

JA'B’ 

=  (1  +  pp‘)~ 

1[—p'Saa  +  pSbb  +  Sab  -  Pp‘Sba)] 

(3.18 

S' 

^B'A' 

=  (l+pp*r 

1  [ — +  pSbb  +  Sba)  —  PP~Sab\ 

(3.19 

S' 

dB'B' 

=  (1 +  PP")- 

l[p'  Saa  +  Sbb  —  p‘(Sab  +  -Sen)] 

(3.20 

The  above  relationships  satisfy  the  transformation  invariants  i.e. 


det  ([S(A,B)])  =  det  ([S’(A'B’)])  =  invariant 


and  finally 

Span  ([S(A,  B)))  =  |5^^|2  +  iSael2  +  |5’s/i|2  +  |Sbs|2  =  p  =  Span([S'(  A'fi')]) 

The  previously  mentioned  CO-POLS  (for  minimum  polarization)  are  obtained 
by  setting  S'A,A,  and/or  S'B,B,  in  equation(3.17  or3.20)  to  zero  and  for  the  monostatic 
case  leads  to  equation(3.21)  shown  below: 


-Sab  +  /  —  J  SAB  —  SbbSaa 
Pi, 2  =  - ? - 

JBB 


(3.21) 


Similarly,  by  setting  S'A,B,  and  Sb'a1  to  zero  the  X-POL  (for  maximum  polarization) 
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is  defined  by  the  following  equation(3.22): 


+  /  —  \Jr\  +  4R2R3 
P  1,2  =  - 


2  R2 


(3.22) 


w 


here 


#1  =  (l-fi’esl2  -  ISUal2  —  Sab^ba  +  SbaS*ab) 

#2  =  +  SbaS'4A) 

R3  =  ( SaaS'ab  +  SabSbb ) 

The  latitude  and  longitude  for  both  the  CO-POL  and  the  X-POL  are  deter¬ 
mined  via  equation(3.23)  and  equation(3.24)  as  shown  below: 


cos(0)  = 


1  -  M2 

1  =  M5 


(3.23) 


and 


tan(<^  )  =  — 


i™\q\ 

Re\q\ 


(3.24) 


The  variable  q  in  equation(3.23  and  3.24)  is  defined  by  the  following  equation(3.25): 


_  1  -jp 
9  1  +3P 


(3.25) 


The  operation  and  other  relevant  facets  of  the  Mathematica  (m  routine  are 
documented  within  the  programme  a  copy  of  which  appears  at  Appendix  B.  Upon 
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initialization,  the  Mathematica  tm  routine  automatically  reads  in  data  files  and  any 
required  operative  programmes,  for  example  parametric  or  animation  routines.  Pro¬ 
gramme  output  consists  of  Poincare  data  and  polarization  plots  which  reflect  the 
polarization  behaviour  of  the  canonical  within  the  resolution  cell.  Examples  of  this 
output  are  discussed  in  the  results  section  of  this  thesis. 

3.3.5  Polanmetric  Data  Within  a  Cell.  The  above  theory  and  associated 
software  provided  a  ready  and  visible  means  whereby  intra  resolution  cell  polari- 
metrics  could  be  examined.  The  actual  methodology  employed  is  contained  in  the 
remainder  of  this  chapter. 

The  first  phase  of  the  experiment  consisted  of  verifying  that  a  pure  two  bounce 
canonical  (dihedral)  and  a  pure  single  bounce  canonical  (flat  plate)  exhibited  the  cor¬ 
rect  relationships  between  their  individual  optimal  polarizations.  This  was  achieved 
by  '‘placing'’  each  of  the  shapes  centrally  in  the  resolution  cell  and  examining  the  re¬ 
sulting  output.  Subsequent  to  this  phase  separate  experiments  were  conducted  with 
each  of  the  canonicals  to  determine  whether  or  not  the  polarization  characteristics 
were  a  function  of  the  shapes’  position  within  the  resolution  cell.  This  involved  mov¬ 
ing  the  canonicals  in  all  three  dimensions  and  then  examining  the  resulting  series  of 
data  for  corelation.  The  actual  paths  that  the  shapes  were  moved  along  are  shown 
in  Chapter  IV. 

The  second  phase  of  the  experiment  entailed  placing  the  two  canonicals  within 
a  single  resolution  cell  and  examining  the  resulting  composite  polarization  charac- 
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teristics.  Initially  the  geometric  center  of  the  two  canonicals  was  placed  in  the  center 
of  the  cell  and  the  returns  were  analyzed  in  an  attempt  to  ascertain  whether  their 
combined  prescence  could  have  been  determined  without  apriori  knowledge.  Subse¬ 
quent  to  this  the  composite  shape  was  moved  around  within  the  cell  to  determine 
any  positional  dependency  of  their  combined  polarization  signature.  The  final  seg¬ 
ment  of  this  phase  involved  simultaneously  moving  each  shape  along  a  different  path 
within  the  resolution  cell  to  ascertain  the  dependence  of  the  polarization  character¬ 
istics  on  canonical  speration.  The  exact  intra  cell  trajectories  used  are  contained  in 
ChapterIV. 

The  final  phase  of  the  experiment  involved  changing  the  relative  sizes  of  the 
two  canonicals  in  an  attempt  to  ascertain  whether  there  was  a  size  threshold  above 
which  one  of  the  canonicals’  polarization  characteristics  would  dominate  the  cell's 
polarization  return.  The  size  of  the  canonicals  used  in  this  and  the  two  previous 
phases  are  detailed  in  ChapterIV. 

In  each  of  the  above  phases  SarTooltm  was  configured  such  that  both  the  az¬ 
imuthal  scan  was  limited  to  ±  45  degrees  either  side  of  the  shapes’  center  while  the 
"elevational  angle”  was  held  constant  at  90  degrees.  Furthermore  in  all  cases  the 
canonicals  were  held  “face  on”  to  the  transmitter  as  illustrated  in  Figure  3.4 
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Figure  3.4.  SarTooF”1  Target  Orientation 


3.4  Perceptron  Analysis 

At  the  conclusion  of  the  above  experiments  an  investigation  into  the  appli¬ 
cability  of  perceptron  networks  to  the  task  of  target  classification  was  conducted. 
This  investigation  was  not  intended  to  be  all  encompassing  but  was  designed  to  pro¬ 
vide  an  initial  insight  into  the  potential  of  coupling  polarimetric  data  to  perceptron 
network  analysis.  Accordingly,  the  next  few  paragraphs  will  briefly  summarize  the 
characteristics  of  the  two  perceptron  networks  used  during  this  thesis  effort. 

3-4-1  Single  Level  Perceptron.  This  network  was  based  on  an  IEEE  article 
written  by  Lippman  [23]  in  which  he  presents  a  perceptron  convergence  procedure 
originally  developed  by  Rosenblatt  [24].  Rosenblatt’s  algorithm  is  based  on  the 
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single  layer  perception.  This  type  of  network  can  employ  continuous  as  well  as 


discrete  inputs  and  is  capable  of  quickly  recognizing  simple  patterns.  The  algorithm 
developed  by  Rosenblatt  involved  initializing  the  connection  weights  and  threshold 
values  to  small  random  non  zero  values.  The  perceptron  is  then  presented  with 
the  classes  of  interest  and  is  prompted  to  make  a  classification.  When  the  network 
returns  an  incorrect  classification  the  weights  are  updated,  otherwise  they  remain 
unchanged.  Accordingly,  over  a  period  of  time,  the  weights  converge  such  that 
the  network  achieves  its  maximum  classification  performance.  A  gain  term  exists 
within  the  algorithm  which  allows  the  operator  to  opt  for  faster  adaptation  weights 
at  the  expense  of  weight  stability.  The  resulting  network  is  ultimately  capable  of 
constructing  a  hyperplane  which  bisects  the  two  separable  classes  thereby  facilitating 
simple  class  classification.  Figure  3.5  shows  the  basic  perceptron  unit  and  associated 
hyperplane.  A  program  was  written  which  incorporated  Rosenblatt’s  algorithm  and 
was  used  as  a  classification  aid  during  the  experimental  analysis  section  of  this  thesis. 
A  copy  of  this  program  and  some  sample  outputs  appears  at  Appendix  B. 

3-4-2  Multi  Layer  Perceptron.  As  mentioned  above  the  single  level  percep¬ 
tron  is  adept  at  recognizing  simple  and  readily  separable  data.  Unfortunately,  the 
network  experiences  much  difficulty  when  the  decision  regions  of  the  classes  overlap. 
Accordingly,  the  more  complex  experimental  data  was  used  as  an  input  to  a  multi 
layer  perceptron  network.  The  network  used  was  developed  by  G.Tarr  [25]  and  em¬ 
ployed  Werbose’s  back  propagation  algorithm  with  a  momentum  term  included,  as 
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described  in  Lippman’s  article  [23].  The  resulting  program  is  an  extremely  useful 
piece  of  software  which  allows  the  operator  to  structure  the  network  according  to 
the  level  and  complexity  of  the  input  data.  The  network  used  in  this  thesis  effort 
consisted  of  56  input  nodes  and  two  hidden  layers.  The  first  and  second  hidden  layer 
consisted  of  20  and  10  nodes  respectively.  Further  details  regarding  this  network 
appear  in  Chapter  4. 

3.5  Summary 

The  methodology  adopted  in  this  thesis  involved  understanding  the  limita¬ 
tions  of  non  polarimetric  SAR  resolution  and  then  designing  a  software  environment 
conducive  to  the  examination  of  polarimetric  SAR  data.  The  examination  of  the 
polarimetric  data  was  ultimately  directed  at  ascertaining  whether  polarimetric  an- 
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alyzes  could  achieve  intra  cell  identification  without  apriori  knowledge.  Perceptron 
networks  were  used  during  this  process.  The  results  of  the  above  endeavours  are 
contained  within  Chapter  4  of  this  report. 
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IV.  Results 


4.1  Introduction 

As  noted  in  Chapter  3,  this  thesis  was  intended  to  investigate  three  primary 
areas  of  concern.  These  areas  are  as  detailed  below: 

1.  The  validation  of  SarTool'”1. 

2.  The  development  of  a  composite  software  package. 

3.  The  investigation  of  intra  resolution  cel!  polarimetrics. 

In  addition  to  the  above  primary  areas  of  interest  the  collected  polarimetric 
data  was  processed  and  then  used  as  input  to  two  perceptron  networks  in  an  at¬ 
tempt  to  gauge  the  potential  of  Artificial  Intelligence  (AI)  in  the  discrimination  of 
polarimetric  patterns.  The  results  of  the  AI  endeavour  and  of  the  primary  areas  of 
investigation  are  described  in  the  following  sections  of  this  chapter. 

4-2  The  Validation  of  SarToof™ 

The  development  of  SarToolfm  is  an  ongoing  process  and  as  such  the  sponsor 
of  this  thesis  needed  to  know  how  well  the  existing  documentation  correlated  with 
the  currently  released  version  of  the  software.  Because  an  exhaustive  validation  of 
the  software  and  documentation  was  beyond  the  scope  and  time  limitations  of  this 
thesis  effort,  it  was  agreed  to  only  investigate  the  dihedral  and  flatplate  canonicals 
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for  correspondence  between  theoretical  results,  documented  algorithms  and  software 
performance.  The  basis  for  the  validation  was  chosen  as  being  Griesser's  article  ti¬ 
tled  Backscatter  Analysis  of  Dihedral  Corner  Reflectors  Using  Physical  Optics  and 
the  Physical  Theory  of  Diffraction  [17]  in  which  he  presents  the  theoretical  and  mea¬ 
sured  data  of  several  dihedrals.  Griesser  used  a  9.4GHZ  radar  on  a  dihedral  with 
dimensions  of  5.6086  A.  The  Equations  shown  in  [26]  were  coded  using  MATHCAD"71 
and  Mathematical  and  were  then  compared  to  those  shown  in  Griesser’s  article. 
In  both  cases  there  were  marked  differences  between  Griesser’s  article  and  the  out¬ 
put  predicted  by  the  SarTooltm  documentation.  Investigations  with  TASC  revealed 
thai  the  original  equations  for  the  dihedral  cited  in  the  TASC  documentation  were 
incorrect  although  the  actual  software  returned  the  correct  prediction. 

A  further  anomaly  with  the  SarTooltm  package  was  discovered  during  the  pre¬ 
liminary  stages  of  experiment  design.  The  version  of  SarTooltm  held  by  the  sponsor 
failed  to  acknowledge  the  presence  of  canonicals  within  a  resolution  cell  unless  the 
canonicals  were  in  excess  of  15  A  (A  being  the  wavelength  of  the  target  radar).  TASC 
acknowledged  this  inadequacy  and  released  a  software  patch  which  apparently  cor¬ 
rects  the  problem. 

One  final  point  of  interest  regarding  the  performance  of  SarToolim  is  that  when 
two  canonicals  are  within  a  single  resolution  cell  the  software  does  not  make  provision 
for  any  interactions  between  the  two  shapes.  This  results  in  the  cell  response  being 
merely  an  addition  of  the  individual  responses  of  the  two  canonicals.  This  revelation 


4-2 


will  be  discuss*-'!  in  greater  detail  in  the  following  section  of  this  chapter. 

The  above  observations  regarding  the  accuracy  or  limitations  of  SarTool(m  are 
not  intended  to  be  extrapolated  into  a  general  judgement  of  the  entire  software 
package.  Nevertheless,  even  this  limited  investigation  did  reveal  several  unexpected 
anomalies  and  certainly  indicates  that  any  future  experiments  or  thesis  efforts  in¬ 
volving  SarTooi(rn  should  be  preceded  by  an  examination  of  the  involved  software 
modules  and  their  documentation.  It  should  also  be  noted  that  TASC  personnel  were 
most  helpful  on  every  occasion  when  these  problems  were  brought  to  their  attention. 

4-3  The  Development  of  a  Composite  Software  Package 

A  requirement  of  this  thesis  was  to  develop  a  software  package  capable  of 
interfacing  with  SarToolJm  such  that  SAR  data  could  be  extracted  and  used  to 
investigate  radar  nolarimetric  behavior  within  a  single  resolution  cell.  Appendix  B 
contains  listings  of  the  software  developed  and  used  during  this  thesis  e.  rt.  The 
following  paragraphs  discuss  the  functions  of  each  of  these  software  modules  and 
illustrates  how  they  interact  to  form  a  composite  data  collection  software  package. 

4-3.1  Data  Collection  Procedure.  Initially  several  directories  were  established 
on  the  SUN  workstation  corresponding  to  each  of  the  major  experiments  described  in 
the  following  section.  Within  each  of  these  directories  were  sub  directories  dedicated 
to  the  various  combinations  of  geometric  shapes  being  examined.  Each  of  the  sub 
directories  contained  the  following  three  files: 
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1.  SarTool(rn  .dat  File. 


2.  Canonical  Reflector  File. 

3.  \lodTfi.dat  file. 

Tne  fast  two  oi  the  above  files  are  integral  parts  of  the  SarTool<m  package  and 
as  such  are  described  in  Appendix  A.  The  third  file  i.e.  the  modrfi.dat  file  defines  the 
start  and  finish  position  of  the  canonicals  and  the  number  of  “steps”  or  increments 
taken  between  those  two  positions.  The  modTji.dat  file  is  written  by  the  operator 
according  to  the  required  experiment  configuration  and  is  used  as  an  input  to  the 
niodr/i  software  described  in  Appendix  A.  SarToolim  is  invoked  at  the  start  position, 
the  stop  position  and  at  every  increment  position  A  sample  rv.odTf<.dat  file  is  shown 
in  Table(4.1). 


Table  4.1.  Example  of  a  ModTfi.dat  File 


'dihed’ 

NAME  OF  INPUT  REFLECTOR  FILE 

TO’ 

NUMBER  OF  INCREMENTS 

’2’ 

CANONICAL  IDENTIFIER 

’X,Y,Z’ 

START  POSITION  OF  CANONICAL 

X,Y,Z’ 

FINISHING  POSITION  OF  CANONICAL 

The  last  three  lines  of  Table(4.1)  can  be  repeated  up  to  100  times  thereby 
allowing  a  maximum  of  100  canonicals  to  be  simultaneously  moved  through  the 
radar’s  field  of  view.  As  mentioned  earlier  one  of  the  aims  of  this  thesis  was  to 
facilitate  the  ready  acquisition  of  polarimetric  data;  accordingly,  the  data  collection 
process  was  automated  as  much  as  possible  and  consists  of  the  following  steps: 
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1.  Step  !.  In  each  SUN  experiment  sub  directory  there  are  initially  three  files 
i.e.  the  SarTool'm.dat  file,  the  Canonical  Reflector  File,  the  Modrfl.ilat  file 
and  a  shell  program  entitled  “experiment. sh”.  At  the  commencement  of  each 
experiment  the  shell  program  is  invoked  which  then  performs  the  following 
functions: 

(a)  Runs  inodTj\  thereby  creating  formatted  reflector  files  corresponding  to 
the  input  parameters  of  the  modrji.dat  file. 

(b)  Converts  the  formatted  reflector  files  to  binary  files. 

(c)  Copies  the  above  binary  files  into  a  temporary  reflector  files  a  Inch  are 
subsequently  “called”  by  the  Sar.dat  file. 

(d)  Invokes  SarTool(m  thereby  resulting  in  a  rpdb  file  for  each  “step"  position 
contained  within  the  modTji.dat  file. 

(e)  Rim.'  r  he  dbstrip  program  which  outputs  stripped  “in  phase”  and  "quadra¬ 
ture''  data  corresponding  to  each  of  the  rpdb  files.  The  series  of  flies  con¬ 
taining  the  “stripped”  data  are  titled  "testout.O"  through  to  “testout.X" 
where  X  is  the  number  of  increments  taken  through  the  resolution  cell. 

(f)  Lra-.es  all  of  the  reflector  files  to  facilitate  the  next  run. 

2.  Step  2.  All  of  the  above  actions  associated  with  the  shell  program  are  “trans¬ 
parent"  to  the  operator.  However,  when  the  shell  program  has  finished  running 
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the  resulting  "test"ut  output  files  must  he  transferred  to  the  Macintosh  com¬ 
puter  via  File  Transfer  Protocol  and  NCSA  Telnet 

3.  Step  3.  Once  the  test  output  files  are  transferred  to  the  Macintosh  thev  are 
individually  processed  by  the  Mathematical  ’•ccjtinc  with  the  resulting  output 
being  stored  within  individual  experiment  sub  directories. 

4.  Step  4 •  The  output  from  the  Mathematical  routine  is  then  converted  into 
an  appropriate  form  for  input  into  either  the  single  level  percept ron  or  into 
G.Tarr’s  multi  level  perceptron  network.  Conversion  for  Tarr's  network  is 
performed  by  the  “convert'1  program  which  is  listed  in  Appendix  B. 

4-3.2  Ancillary  Software.  The  above  procedures  required  both  coupling  be¬ 
tween  the  SUN  workstation  and  the  Macintosh  computer  and  the  use  of  a  transfer 
protocol,  for  example,  NCSA  Telnet.  Moreover,  an  IBM  (AT)  Personal  Computer 
was  used  quite  extensively  as  a  development  tool,  thereby  necessitating  the  use  of 
the  TOPS  Translator.  Input  file  editing  was  done  via  a  combination  of  Macintosh 
generic  software  and  the  VI  editor. 

A  diagrammatic  representation  showing  the  integration  of  the  composite  soft¬ 
ware  package  is  shown  in  Figure  4.3.2. 

4-4  Intra-Cell  P olarimetrics 

The  investigation  of  Intra-Cell  polarimetrics  was  centered  around  four  basic 
experimental  blocks.  A  discussion  of  these  experimental  blocks  and  an  interpretation 
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Figure  4.1.  Integration  of  Software  Packages 


of  their  results  is  contained  in  the  following  paragraphs. 

4-4- 1  Experimental  Block  One.  Experimental  block  one  consisted  of  three 
sub  experiments.  The  objective  of  these  sub  experiments  was  to  ascertain  whether 
or  not  SarTooP771  made  any  distinction  based  on  the  position  of  the  canonical  within 
the  resolution  cell.  To  determine  this  the  SarTooPm.dat  file  was  configured  as  shown 
in  Table  4.4.1. 

In  addition  to  the  above  configuration  the  resolution  cell  was  given  a  nomi¬ 
nal  height  of  2m  thereby  resuming  in  a  2m  cube  resolution  cell.  This  cell  and  the 
numbering  convention  adopted  are  shown  in  Figure  4.4.1. 

Sub  experiment  one  was  divided  into  three  sections  as  described  overleaf: 
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Table  4.2.  SarTool(m  Sensor  and  System  Data  (Sartool.dat) 


'expl' 


'dihed' 


'O’ 


'O' 


'O' 


•9.4  GHZ’ 


’2.0’ 


’2.0' 


’.01’ 


TITLE  OF  SIMULATION  RUN 


SPECIFIES  INPUT  REFLECTOR  FILE 


X  POSITION  OF  SENSOR 


Y  POSITION  OF  SENSOR 


Z  POSITION  OF  SENSOR 


FLIGHT  PATH  ANGLE  OF  SENSOR 


DIVE  ANGLE  OF  SENSOR 


AVERAGE  SENSOR  SPEED  (SAR  MODE) 


X  POSITION  OF  TARGET  REL  TO  MAP  CENTER 


Y  POSITION  OF  TARGET  REL  TO  MAP  CENTER 


Z  POSITION  OF  TARGET  REL  TO  MAP  CENTER 


EULER  ANGLE  PHI  FOR  TGT  FRAME  ROTN 


EULER  ANGLE  THETA  FOR  TGT  FRAME  ROTN 


EULER  ANGLE  PSI  FOR  TGT  FRAME  ROTN 


TRANSMITTED  FREQUENCY 


ELEVATN  3db  BEAMWIDTH  (DEG) 


AZIMUTH  3db  BEAMWIDTH  (DEG) 


ELEVATION  BORESIGHT  ANGLE  (DEG) 


’90.0’ 

AZIMUTH  BORESIGHT  ANGLE  (DEG) 

’2.0’ 

RANGE  RESOLUTION  (m) 

’2.0’ 

AZIMUTH  RESOLUTION  (m) 

T 

NUMBER  OF  RANGE  BINS 

T 

NUMBER  OF  AZIMUTH  BINS 

’200’ 

NUMBER  OF  PULSES  TRANSMITTED 

T 

TARGET  DATABASE  RECALCULATION  I/VAL 

T 

TERRAIN  GENERATION  OPTION  (SAR) 

'O' 

AVERAGE  TERRAIN  RCS  (dbsm) 

Y 

SAR  SPOLTLIGHT  OPTION 

V 

ISAR  MODE  OPERATION 

’O’ 

ASPECT  ANGLE  INCREMENT  FOR  QKTOOL 

T 

MOTION  COMPENSATION  OPTION 

T 

REFLECTOR  THINNING  OPTION 

'O' 

THINNING  THRESHOLD  VALUE  (dbsm) 

Y 

TARGET  SHADOWING  OPTION 

Y 

TERRAIN  SHADOWING  OPTION 

T 

REFLECTOR  I/P  REPORT  OPTION 

V 

REFLECTOR  D/BASE  REPORT  OPTION 

Y 

RANGE  PROFILE  OUTPUT  OPTION 
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Figure  4.2.  Experimental  SAR  Resolution  Cell 


1.  Dihedral.  A  dihedral  having  sides  of  length  8.944  cm  was  moved  from  corner 
1  to  corner  6  in  10  steps. 

2.  Flatplate.  A  square  flatplate  (17.88cm)  was  moved  in  10  increments  along  the 
same  path  as  the  dihedral  above. 

3.  Combination.  The  two  canonicals  mentioned  above  were  combined  and  moved 
together  in  10  steps  from  corner  1  to  corner  6. 

Sub  experiments  two  and  three  used  the  same  sized  canonicals  and  were  divided 
into  the  same  three  sections  as  shown  above.  However,  in  sub  experiment  2  the  path 
taken  was  from  corner  eight  to  corner  three  while  in  sub  experiment  3  the  path 
selected  was  from  A  to  B.  Samples  of  the  relevant  Poincare  sphere,  co-polarization 
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and  span  plots  taken  during  Experiment  block  1  are  shown  in  Figures  4.3  to  4.11. 
Note  that  the  figure  nomenclature  “x/y”  signifies  experiment  x,  sub  experiment  y. 
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•  C/Q 


0  *  1  1.5  2.5 


Latitude 

Figure  4.4.  Exp  1/1  -  Dihedral  (Path  1  to  6)  Poincare  Plot 


0  5  10  15  20  25  30 

Pulse  Number 

Figure  4.5.  Exp  1/1  -  Dihedral  (Path  1  to  6)  Span  Plot 
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Figure  4.6.  Exp  1/2  -  Flatplate  (Path  8  to  3)  Poincare  Sphere 
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Figure  4.7.  Exp  1/2  -  Flatplate  (Path  8  to  3)  Poincare  Plot 
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Fi<rniv>  4  8  Exp  1/2  -  Flat  plate  (Path  8  to  3)  Span  Plot 


bO  • 


Figure  4.10.  Exp  1/3  -  Combination  (Path  A  to  B)  Poincare  Plot 

25  i  i  i  i  i  i  I 


20 


0  5  10  15  20  25  30 

Pulse  Number 

Figure  4.11.  Exp  1/3  -  Combination  (Path  A  to  B)  Span  Plot 
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Examination  of  the  of  the  outputs  from  experiment  one  showed  that  different 
canonicals  resulted  in  clearly  different  Poincare  sphere  plots.  Moreover,  the  cell's 


response  was  found  to  be  constant,  regardless  of  the  canonical’s  position  within  the 
cell. 


4-4-2  Experimental  Block  Two.  Experimental  block  two  consisted  of  three 
sub  experiments  which  were  designed  to  measure  the  interaction  of  one  canonical  on 
another  canonical  when  both  the  canonicals  are  co-located  within  the  same  resolution 
cell.  The  three  sub  experiments  are  described  below: 

1.  Exp2  Sub  Experiment  1.  In  this  experiment  a  dihedral  was  moved  from  corner 
2  to  corner  6  in  five  increments  while  at  the  same  time  a  flatplate  was  moved 
from  corner  8  to  corner  4  in  the  same  number  of  increments. 

2.  Exp2  Sun  Expeumcnt  2.  This  sub  experiment  was  a  virtual  repeat  of  the  above 
sub  experiment  except  that  the  canonical  paths  were  reversed  i.e.  the  dihedrad 
was  moved  from  corner  8  to  corner  4  while  the  flatplate  was  moved  from  corner 
2  to  corner  6. 

3.  Exp2  Sub  Experiment  3.  In  this  experiment  one  dihedral  was  moved  from 
corner  2  to  corner  6  while  another  dihedral  was  moved  form  corner  8  to  corner 

4. 

Sample  outputs  from  Experiment  2  are  shown  in  Figures  4.12  to  4.20. 
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Figure  4.12.  Exp.  2/1  Dihedral  Path  (2  to  6)  Flatplate  Path  (8  to  3)  Poincare 
Sphere 


Latitude 


Figure  4.13.  Exp.  2/1  Dihedral  (Path  2  to  6)  Flatplate  (Path  8  to  3)  Poincare  Plot 


Figure  4.15.  Exp.  2/2  Dihedral  (Path  8  to  4)  Flatplate  Path(2  to  6)  Poincare 
Sphere 
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Figure  4.17.  Exp.  2/2  Dihedral  (Path  8  to  4)  Flatplate  (Path  2  to  6)  Span  Plot 
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igure  4.19.  Exp.  2/3  -  Dihedral  (Path  2  to  6)  Dihedral(Path  8  to  4)  Poincare  Plot 
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Figure  4.20.  Exp.  2/3  -  Dihedral  (Path  2  to  6)  Dihedral(Path  8  to  4)  Span  Plot 

Analysis  of  the  output  from  experiment  2  revealed  that  the  SarToo)tTn  algo¬ 
rithms  do  not  make  any  provision  for  interaction  between  the  two  canonicals.  This 
was  evident  in  that  the  cell’s  response  at  the  start  positions  of  each  sequence  was 
identical  to  the  cell’s  response  at  each  subsequent  interval.  This  is  an  important 
revelation,  the  impact  of  which  will  be  discussed  in  a  later  section  of  this  thesis. 
Also  of  interest  in  experiment  block  2  were  the  Poincare  patterns  which  resulted 
from  combining  two  canonicals  in  a  single  resolution  cell.  These  patterns  showed 
clearly  discernible  traits  of  each  primitive’s  co-polarization  signature.  This  result 
influenced  the  conduct  of  future  experiments  and  is  further  detailed  in  the  following 
paragraphs. 
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Table  4.3.  Experiment  Block  Three:  Dihedral  and  Flatplate  Sizes 


1  Dihedral  Size  (cm) 

Flatplate  Size  (cm) 

Exp  No. 

35.776 

37.888 

3a 

53.664 

37.888 

3  b 

71.552 

37.888 

3c 

89.440 

37.888 

3d 

8.944 

37.888 

3e 

5.963 

37.888 

3f 

4.472 

37.888 

3g 

3.578 

37.88 

3h 

4-4-3  Experimental  Block  Three.  This  section  of  the  experiment  was  devoted 
to  determining  the  effect  of  changing  relative  canonical  sizes  on  the  Poincare  sphere 
output.  The  experimental  procedure  involved  holding  the  size  of  the  flatplate  con¬ 
stant  whilst  first  increasing  and  then  decreasing  the  size  of  the  dihedral.  The  actual 
sizes  used  during  this  experiment  are  detailed  in  Table  4.4.3. 

The  outputs  from  experiments  3a,  3e  and  3h  are  shown  in  Figures  4.21  to  4.29 
respectively. 
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Figure  4.21.  Exp.  3a  Combination  -  Poincare  Sphere 
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Figure  4.22.  Exp.  3a  Combination  -  Poincare  Plot 
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Figure  4.25.  Exp.  3e  Combination  -  Poincare  Plot 
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Figure  4.26.  Exp.  3e  Combination  -  Span  Plot 


4-24 


Figure  4.27.  Experiment  3h  -  Poincare  Sphere 
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Figure  4.28.  Exp.  3h  Combination  -  Poincare  Plot 
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Figure  4.29.  Exp.  3h  Combination  -  Span  Plot 


The  above  results  proved  to  be  of  significant  interest.  The  pattern  in  experi¬ 
ment  3a  was  consistent  with  previous  results  in  that  it  clearly  showed  the  dominance 
of  the  dihedral  co-polarization  nulls  over  those  of  the  flat  plate.  However,  in  ex¬ 
periments  3e  and  3h  the  flatplate  is  assuming  dominance  due  to  its  increased  size 
over  the  dihedral  (see  Table  4.4.3).  Accordingly,  there  is  a  marked  movement  of 
the  co  polarization  nulls  towards  the  “poles”  of  the  sphere.  A  point  of  particular 
interest  revealed  by  experiment  block  three  was  that  this  pattern,  although  different 
tor  various  size  combinations,  appeared  to  have  a  consistent,  predictable  trend.  This 
trend  is  evidenced  by  noting  the  similarities  between  the  Poincare  sphere  results  of 
experiment  3e  and  experiment  3h.  The  span  of  the  co-polarization  nulls  was  also 
of  interest  in  that  it  showed  a  clear  distinction  between  when  the  dihedral  was  the 
dominant  factor  and  when  the  flat  plate  assumed  dominance.  The  results  of  exper¬ 
iment  3  prompted  an  investigation  into  the  possible  role  that  perceptron  networks 
could  play  in  performing  identification  of  cell  contents  via  an  examination  of  the 
co-polarization  null  patterns.  Accordingly,  in  order  to  investigate  this  possibility 
further  one  final  batch  of  experiments,  Experimental  block  four,  were  performed. 
These  experiments,  numbered  6  through  to  48  and  70  through  105  are  described  in 
the  following  section. 

4-4-4  Experimental  Block  Four.  The  primary  purpose  of  this  block  of  experi¬ 
ments  was  to  collect  sufficient  data  on  dihedrals  and  the  combination  of  dihedrals  and 
flatplates  so  that  tests  could  be  conducted  using  perceptron  networks.  Accordingly, 
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experiments  6  through  49  involved  varying  size  dihedrals  centered  within  the  resolu¬ 
tion  cell  while  experiments  70  through  105  contained  various  sized  combinations  of 
flat  plates  and  dihedrals.  The  exact  configuration  of  these  experiments  is  shown  in 
Table  4.4.4.  Note  that  the  dihedral  size  started  at  5cm  square  in  experiment  6  and 
was  incremented  by  5cm  for  each  subsequent  experiment  such  that  in  experiment 
45  the  dihedral  was  200  cm  square.  The  sizes  of  the  dihedrals  in  experiments  46 
through  48  are  as  shown. 

J,.5  Discussion  of  Poincare  Plots 

Of  note  in  each  of  the  above  experiment  blocks  was  the  "telling”  character¬ 
istics  of  the  Poincare  sphere  plots.  Although  theory  predicts  distinctive  plots  from 
isolated  canonicals  the  level  of  information  retrievable  from  the  sphere  plots,  when 
combinations  of  canonicals  were  involved,  was  not  expected.  This  latent  potential 
for  information  extraction  is  evident  in  most  of  the  discussed  plots  i.e.  even  when 
one  canonical  was  physically  dominant  it  was  unable  to  completly  mask  the  polari- 
metric  effect  of  the  other  canonical.  Specific  examples  of  this  effect  can  be  seen  by 
examining  the  Poincare  plots  of  experiment  3e  and  3h.  In  both  these  experiments 
the  flat  plate  is  physically  the  dominant  canonical  and  yet  the  Poincare  plots  show 
co-polarization  nulls  which  clearly  indicate  the  prescence  of  other  than  a  single  flat 
plate.  Similarly  in  incidences  where  the  dihedral  is  the  physically  dominant  canoni¬ 
cal.  for  example  in  experiment  86,  the  distribution  of  its  co-polarization  nulls  within 
the  Poincare  plot  suggest  that  another  canonical  may  be  present.  In  addition  to 


Table  4.4.  Experiment  Block  Four:  Dihedrals  and  Combinations 


Dihedral  Size  (cm) 


Flatplate  Size  (cm) 


Exp  No. 


5  -  200 
27 
33 
64 
20 
20 
20 
20 
20 
40 
40 
40 
40 
40 
60 
60 
60 
60 
60 
80 
80 
80 
80 
80 
100 
100 
100 
100 
120 
120 
120 
120 
120 
120 
160 
160 
160 
160 
160 


Not  Applicable 


6  -  45 


the  examples  shown  in  this  chapter  there  are  several  experiment  outputs  contained 
in  Appendix  C  all  of  which  show  the  same  type  of  effects  discussed  above.  These 
somewhat  encouraging  observations  led  to  a  preliminary  analysis  of  the  Poincare 
plots  using  perceptron  networks. 

4.6  Perceptron  Network  Analysis 

The  apparent  trend  of  the  above  results,  in  addition  to  those  presented  in  Ap¬ 
pendix  C,  led  to  a  series  of  experiments  regarding  the  applicability  of  perceptron 
networks  to  polarimetric  recognition  tasks.  These  investigations  were  not  intended 
to  be  exhaustive  but  were  designed  to  yield  an  initial  assessment  of  the  viability  of 
employing  A. I.  in  polarimetric  SAR  radars.  The  perceptron  experiments  were  di¬ 
vided  into  two  primary  blocks:  those  involving  a  single  layer  perceptron  network  and 
those  involving  a  multi-level  perceptron  network.  The  results  of  these  experiments 
are  contained  in  the  following  paragraphs. 

4.6.1  Single  Level  Perceptron  Network.  The  Mathematica*"1  output  was  for¬ 
matted  as  latitude  and  longitude  vectors  and  then  used  as  input  to  the  single  level 
perceptron  network  which  appears  in  Appendix  B.  The  network  was  written  such 
that  it  could  either  train  and  test  on  the  same  set  of  data  or  train  on  one  set  of  data 
and  test  on  another  set  of  data.  A  random  number  generator  (tied  to  the  computer 
clock)  was  used  to  randomize  the  initial  perceptron  weights  and  the  selection  of  train 
and  test  vectors.  The  results  gained  using  this  network  are  shown  in  Table  4.6.1. 
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Table  4.5.  Single  Level  Perception  Classification  Results 


Experiment  Number 

Vectors  Used 

Percentage  Right 

1 

ldlf 

100 

2 

2d2f 

100 

3 

3d3c 

52 

4 

ld3c 

48 

5 

ld4c 

54 

Note  that  the  vector  nomenclature  xkxk  decodes  as  shown  below: 


1.  x.  x  equals  the  experiment  block  number  as  discussed  above. 

2.  k.  k  signifies  the  type  of  canonical  responsible  for  the  input  vector:  f  =  flat- 
plate,  d  =  dihedral  and  c  =  combination  flatplate  and  dihedral. 


The  results  of  experiments  1  and  2  above  make  intuitive  sense  in  that  the 
data  representing  a  dihedral  and  a  flatplate  is  clearly  separable  (refer  to  example 
outputs  in  previous  section  of  this  chapter).  Accordingly,  the  perceptron  needed 
very  little  exposure,  for  example  10  vectors  of  each  class,  before  it  could  attain  a 
100%  classification  accuracy. 


The  results  of  experiments  3,4  and  5  above  emphasize  the  degree  of  overlap 
which  exists  between  a  dihedral  and  a  flatplate/dihedral  combination.  This  overlap 
is  evident  in  the  output  examples  contained  in  this  chapter.  The  network  in  experi¬ 
ments  3,4  and  5  was  trained  for  2000  iterations  on  an  equal  number  of  dihedral  and 
combination  vectors  after  which  it  was  tested  on  500  vectors.  The  poor  performance 
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was  not  unexpected  and  indicates  the  inability  of  a  single  layer  network  to  separate 
"single  hit”  complex  data. 

The  above  experiments  suggest  that,  notwithstanding  the  obvious  limitations 
of  single  layer  perceptrons,  they  may  be  suited  to  performing  the  initial  evaluation 
of  cell  contents.  Furthermore  a  single  layer  network  is  clearly  capable  of  providing  a 
"count”  of  single  and  double  bounce  returns  from  a  particular  cell  of  interest.  These 
counts  could  then  be  used  to  estimate  the  contents  of  tlm  cell. 

Although  the  above  single  level  “count”  procedure  may  have  applicability  in  the 
early  stages  of  decision  making  it  does  not  harness  the  sequence  of  patterns  evident 
in  the  Poincare  sphere  plots.  These  patterns  can  only  be  harnessed  by  treating  the 
whole  sequence,  not  just  individual  returns,  as  an  input  to  a  perceptron  network. 
Accordingly,  the  additional  experiments  of  block  four  were  run  and  used  as  inputs 
to  G.Tarr’s  multi-level  perceptron  network. 

4-6.2  Multi  Level  Perceptron  Network.  The  output  of  all  the  dihedrals  and 
dihedral/flatplate  combinations  were  converted  into  a  suitable  form  for  input  into 
Tarr’s  network  by  the  “convert”  program  shown  in  Appendix  B  .  A  total  of  93 
vectors  were  collected,  46  of  these  vectors  were  used  to  train  the  network  while  the 
remaining  47  were  used  as  previously  “unseen”  testing  vectors.  The  actual  network 
configuration  used  consisted  of  56  input  nodes,  two  output  nodes  and  two  hidden 
layers.  The  first  hidden  layer  had  20  nodes  whilst  the  second  layer  consisted  of 
10  nodes.  A  total  of  30  separate  runs  were  made  with  the  initial  weights  being 
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randomly  initialized  and  the  order  of  vector  presentation  being  different  for  each 
individual  network  run.  The  classification  results  of  these  runs  are  shown  in  Figures 
1.6.2  and  4.6.2.  The  distinction  between  “good”  and  “right”  is  that  in  the  former 
case  the  conditional  probability  of  correctly  classifying  the  cell  contents  is  between  .5 
and  .9  while  to  be  gauged  as  “right”  the  conditional  probability  must  be  equal  to  or 
greater  than  .9.  Figure  4.6.2  shows  the  decrease  in  recognition  error  as  the  number  of 
training  iterations  is  increased.  Table  4.6  shows  the  conditional  probabilities  each  of 
the  decision  outputs.  As  in  the  previous  discussion  the  letter  d  represents  a  dihedral 
while  the  letter  c  represents  the  combination  of  a  dihedral  and  a  flatplate. 

The  results  for  the  multi  layer  perceptron  network  are  of  particular  interest 
due  to  their  sharp  training  rise  time  and  the  associated  classification  performance. 
This  indicates  that  the  network  finds  the  two  patterns  readily  separable.  Of  note 
is  that  the  above  performance  was  made  only  on  polarimetric  data  i.e.  consistent 
with  the  scope  of  this  project  the  span  data  was  not  used  as  a  classification  aid. 
The  inclusion  of  span  data  would  almost  certainly  enhance  the  net’s  classification 
performance.  Also,  on  a  cautionary  note,  it  should  be  emphasized  that  all  of  the 
experiments  conducted  during  this  thesis  were  conducted  in  a  “noise  free”  environ¬ 
ment.  Furthermore,  as  mentioned  earlier,  SarTooltm  does  not  cater  for  inter  canonical 
interference. 
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Iterations  x  1000 


Figure  4.31.  Percent  Classified  “Right” 
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Table  4.6.  Conditional  Probability  Table 


PDF(C/C) 

PDF(D/C) 

PDF(C/D) 

P(D/D) 

P(good ) 

netl 

0.958 

0.042 

0.000 

1.000 

0.978 

net2 

0.917 

0.083 

0.000 

1.000 

0.957  1 

net3 

0.917 

0.083 

0.000 

1.000 

0.957 

net4 

0.875 

0.125 

0.000 

1.000 

0.935  ; 

net5 

0.917 

0.083 

0.000 

1.000 

0.957  1 

net6 

0.875 

0.125 

0.000 

1.000 

0.935  | 

net? 

0.958 

0.012 

0.000 

1.000 

0.978 

ne  .8 

0.917 

0.083 

0.000 

1.000 

0.957  ! 

net9 

0.917 

0.083 

0.000 

1.000 

0.957 

net  10 

0.958 

0.042 

0.000 

1.000 

0.978 

netl  1 

0.917 

0.083 

0.000 

1.000 

0.957 

net  12 

0.917 

0.083 

0.000 

1.000 

0.957 

net  13 

0.917 

0.083 

0.000 

1.000 

0.957 

netl4 

0.875 

0.125 

0.000 

1.000 

0.935 

net  15 

0.917 

0.083 

0.000 

1.000 

0.957 

netl6 

0.917 

0.083 

0.000 

1.000 

0.957 

net  17 

0.917 

0.083 

0.000 

1.000 

0.957 

netl8 

0.875 

0.125 

0.000 

1.000 

0.935 

net  19 

0.917 

0.083 

0.000 

1.000 

0.957 

net20 

0.875 

0.125 

0.000 

1.000 

0.935 

net21 

0.875 

0.125 

0.000 

1.000 

0.935 

net22 

0.875 

0.125 

0.000 

1.000 

0.935 

net23 

0.875 

0.125 

0.000 

1.000 

0.935 

net24 

0.917 

0.083 

0.000 

1.000 

0.957 

net25 

0.917 

0.083 

9.000 

1.000 

0.957 

net,26 

0.917 

0.083 

0.000 

1.000 

0.957 

net27 

0.958 

0.042 

0.000 

1.000 

0.978 

net28 

0.917 

0.083 

0.000 

1.000 

0.957 

net29 

0.958 

0  042 

0.000 

1.000 

0.978 

net30 

0.875 

0.125 

0.000 

1.000 

0.935 

Average 

0.911 

0.089 

0.954 

Variance 

0.001 

0.001 

NA 

NA 

0.000 

St  Devr. 

0.028 

0.028 

NA 

NA 

_ 

0  015 

4.4  7 


',.7  Summary 

Notwithstanding  the  provisos  of  the  above  paragraph,  the  results  discussed 
during  this  chapter  provide  some  useful  information  on  the  polarimetric  behavior 
within  a  resolution  cell.  Moreover,  although  SarTool4"1  has  certain  limitations,  it 
can  (when  combined  with  other  software)  provide  a  reasonable  medium  for  the  in¬ 
vestigation  of  radar  polarimetrics.  The  results  gained  during  this  thesis  support  this 
contention  and  seem  to  indicate  that  classification  of  intra  cell  contents  is  possible 
using  only  polarimetric  data.  Finally,  the  use  of  single  and  multi  layer  percept ron 
networks  appears  to  have  much  potential  in  this  classification  process. 
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V.  Conclusions 


5.1  Introduction 

This  thesis  effort  was  directed  at  an  examination  of  the  polarimetric  behavior 
of  canonicals  within  a  resolution  cell.  An  integral  and  necessary  part  of  this  venture 
involved  a  limited  evaluation  of  SarToolfm  and  the  development  of  several  software 
modules  such  that  SarTool<m  could  be  readily  employed  in  future  polarimetric  stud¬ 
ies.  The  conclusions  reached  during  the  course  of  these  endeavors  and  suggestions 
for  future  research  efforts  are  discussed  in  the  following  paragraphs. 

5.2  Evaluation  of  SarToofm 

Time  constraints  and  the  primary  objectives  of  this  thesis  precluded  an  exten¬ 
sive  evaluation  of  SarTooltrn.  Nevertheless,  an  early  investigation  of  SarTooi"  ',  and 
subsequent  discoveries  made  during  the  course  of  this  thesis  allow  the  following  main 
conclusions  to  be  made: 

1.  Documentation.  The  level  and  content  of  SarTooff™  documentation  appears 
to  be  lagging  the  development  of  the  software.  An  example  of  this  observation 
is  the  apparent  confusion  which  exists  over  the  documented,  versus  coded, 
reflection  algorithms. 

2.  Configuration.  Early  problems  were  evident  in  SarTool  Vm  “acknowledgement" 
of  small  canonicals  within  resolution  cells.  Although  this  problem  has  appar- 


ently  been  remedied  there  seemed  to  be  no  structured  way  in  which  the  user 
of  the  code  could  access  certain  software  configuration  details,  for  example 
Version  listings  and  ‘"patch”  status. 

3.  Performance.  The  implementation  of  SarTool(m  is  such  that,  except  for  a 
shadowing  option,  the  algorithms  do  not  cater  for  the  interaction  of  canonicals 
which  are  co-located  within  the  same  resolution  cell. 

The  above  conclusions  are  not  intended  as  criticisms  of  SarToofi”1.  However, 
they  are  included  for  their  potential  worth  in  future  preliminary  stages  of  experiment 
design. 

5.3  Developed  Software 

The  software  developed  during  this  thesis  has  been  exhaustively  described 
within  the  body  of  this  thesis  and  software  listings  appear  in  Appendix  B.  Although 
the  mentioned  ancillary  software  greatly  enhanced  the  use  of  SarTooltm  there  is  po¬ 
tential  for  further  improvements  in  the  automation  of  result  gathering.  At  present 
each  SarTooP"1  run  (i.e.  each  movement  of  a  canonical)  must  be  individually  pro¬ 
cessed  by  Mathematica!m  before  being  converted  for  use  by  a  perceptron  network 
(or  other  post  processing  device).  A  future  worthwhile  endeavour  would  be  to  batch 
process  both  the  inputs  and  the  outputs  to  Mathernatica(m  such  that  each  run  of 
SarTooPm  would  automatically  produce  a  decision  output. 
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5.-1  Recognition  of  Intra  -  Cell  Canonicals 

The  results  gained  during  this  thesis  indicate  that  worthwhile  information  can 
be  extracted  by  examination  of  the  polarization  contents  of  the  Poincare  sphere.  A 
particularly  promising  observation  was  that  the  presence  of  a  dihedral  reflector  can 
be  detected  by  visual  inspection  of  the  Poincare  sphere  even  when  the  dihedral  is 
co-located  with  a  flatplate  of  much  larger  dimensions.  As  the  first  in  a  series  of 
studies  on  this  topic  this  thesis  was  constrained  to  only  examining  the  classic  suigic 
and  double  bounce  canonicals  i.e.  the  flatplate  and  the  dihedral.  Notwithstanding 
this  constraint  the  software  that  was  developed  throughout  this  thesis  can  be  applied 
equally  well  to  the  investigation  of  any  canonical  within  single  or  multiple  resolution 
cells.  Moreover,  the  investigation  of  canonicals  where  S^gandSBA  are  not  always 
equal  to  zero  would  facilitate  incorporating  the  cross  polarization  nulls  as  well  as  the 
co-polarization  nulls  into  the  identification  process. 

A  logical  extension  to  the  Poincare  sphere  part  of  this  thesis  would  be  to 
ultimately  build  targets  based  on  as  many  as  necessary  of  the  13  canonicals  modelled 
within  SarToolim.  The  resulting  outputs  could  then  be  processed  by  some  form  of  A. I. 
(perceptron  networks  are  discussed  in  the  next  section)  in  an  attempt  to  gain  target 
recognition.  As  indicated  during  the  previous  chapter  this  thesis  was  conducted 
in  a  ‘‘noise  free”  environment  and  within  certain  practical  limitations  of  SarTool'm; 
accordingly,  any  extension  to  this  effort  should  include  validation  of  the  theoretically 
predicteu  result  within  an  anechoic  chamber. 


5.5  Perceptron  Networks  in  Polarimetrics 

The  use  of  single  and  multi  layer  perceptron  networks  in  the  latter  stages  of  this 
thesis  effort  proved  to  be  of  much  interest.  Accordingly,  the  following  paragraphs  will 
briefly  summarize  the  findings  made  regarding  these  networks  and  will  suggest  areas 
of  future  potential  coupling  between  perceptron  networks  and  polarimetric  studies. 

5.5.1  Single  Level  Perceptron.  The  single  perceptron  network  required  little 
training  before  it  was  able  to  perform  extremely  well  in  the  identification  of  single 
canonicals.  However,  even  with  extensive  training,  the  same  network  was  clearly 
incapable  of  correctly  identifying  a  combination  of  canonicals.  This  inability  arises 
from  the  methodology  used  by  the  single  perceptron  network  to  process  the  cell 
returns.  The  single  perceptron  examines  each  individual  return  from  the  cell  and 
classifies  that  return  as  originating  from  either  a  single  or  double  bounce  canonical. 
This  process  eliminates  the  possibility  of  any  true  pattern  recognition  being  con¬ 
ducted  on  the  returns.  Nevertheless,  useful  information  regarding  ceil  contents  may 
be  extracted  by  even  a  single  layer  perceptron  based  on  relative  counts  of  single  and 
double  bounce  returns.  The  network  contained  in  Appendix  B  would  be  ideal  for 
this  application  as  it  can  maintain  an  ongoing  count  of  the  processed  data. 

5.5.2  Multi  Layer  Perceptron.  The  preliminary  results  gained  using  a  multi 
level  perceptron  network  for  decision  making  appear  to  be  extremely  promising. 
With  little  training  the  multi  level  network  was  able  to  gain  very  high  levels  of 


correct  identification  (approximately  97%)  using  only  co-polarization  polarimetric 
data.  Should  cross  polarization  and  span  data  be  incorporated  into  the  decision 
process  it  seems  possible  that  worthwhile  decision  analysis  could  be  obtained  even 
if  the  library  of  choices  was  increased  to  cover  multiple  combinations  of  canonicals. 

Ultimately  a  structured  approach  involving  a  combination  of  single  and  multi 
layer  networks  may  prove  to  be  computationally  efficient  way  to  perform  analyzes 
of  cell  content.  Initial  decisions  regarding  the  presence  of  single  canonicals  could 
be  made  very  quickly  by  single  layer  networks  while  the  multi  layered  networks 
could  be  used  to  process  the  more  complex  returns  characteristic  of  combinations  of 
canonicals.  Implicit  in  the  above  discussion  is  that  SAR  processing  for  this  type  of 
analyzes  would  need  to  be  based  on  “multiple  looks”  at  each  resolution  cell  rather 
than  on  the  traditional  “one  hit”  approach. 

5.6  Summary 

The  aims  of  this  thesis  effort  were  met:  the  potential  SarTool!m  package  to 
examine  polarimetric  data  was  harnessed  by  ancillary  software  and  it  ’tiple  experi¬ 
ments  examining  cell  polarimetrics  were  conducted.  In  addition  an  initial  investiga¬ 
tion  on  the  applicability  cf  perceptron  networks  to  the  polarimetric  decision  process 
was  conducted.  Logical  extensions  to  this  thesis  include  a  theoretical  and  measured 
study  of  multiple  combinations  of  canonicals.  Finally,  the  incorporation  of  all  po- 
larimetric  and  span  data  into  a  specifically  tailored  perceptron  network  would  seem 


to  be  an  extremely  promising  avenue  of  further  research. 
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Appendix  A.  Relevant  Theory  of  SarTool tm 


A.l  Introduction 

SarTool*”1  is  a  software  package  which  simulates  the  behavior  of  a  polarimefric 
(  HH,  HV,  VH,  YY )  SAR  thereby  allowing  the  scattering  behavior  of  several  canonical 
shapes  to  be  examined.  This  Appendix  provides  a  brief  overview  of  SarTool*”1  as  it 
applies  to  this  thesis  effort. 

A.  2  Background 

SarTool*  m  was  developed  for  the  USAF  by  The  Analytic  Sciences  Corporation 
(TASC).  The  original  contract  number  was  F33615-87-C-  1404  while  the  performance 
of  the  system  is  described  in  a  TASC  Technical  Information  Memorandum  (TIM) 
numbered  5417-  '_’}}( [26]).  The  reader  should  be  aware  that  SarTool*”1  is  still  sub¬ 
ject  to  ongoing  development;  accordingly,  updated  versions  of  the  software  and  the 
associated  documentation  are  periodically  released. 

A. 3  Overview  of  SarTooF™  Pveration 

SarTool*”1  is  intended  to  provide  Electromagnetic  Prediction  of  target  signa¬ 
tures.  As  such  SarTool*"1  uses  Physical  Optics  modelling  to  predict  the  interaction 
between  certain  canonical  geometric  shapes  and  a  fully  polarimetric  SAR.  The  "li¬ 
brary”  of  geometric  shapes  consists  of  thirteen  canonicals  including  elliptical  base 
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t « ■ ; >  hat  reflectors,  cone  reflectors,  ogive  reflectors,  cavity  reflectors,  complex  cavity 
reflect ors.  elliptical  plate  reflectors,  polygonal  base  top  hat  reflectors,  cylinder  re- 
i!'"  tors,  curved  dihedral  reflectors,  ellipsoid  reflectors,  trihedral  reflectors,  dihedral 
r> ■Hectors  and  polygonal  plates.  The  use  of  Solid  Objects  Modelling  allows  these 
■Tapes  to  represent  most  targets  of  interest.  The  "normal"  output  of  SarTooltm  is  a 
SAR  image  the  shape  and  characteristics  of  which  depend  on  the  shapes  illuminated 
i>y  the  radar. 

1  he  user  "interacts"  with  the  SarTool(m  software  via  a  number  of  input  files 
which  are  described  in  the  following  paragraphs. 

A. 3.1  SarTooFm .dat  file.  The  SarTool(m.dat  file  is  the  "master"  input  file 
which  essentially  configures  the  relationship  between  the  radar  and  the  targets.  In- 
pic  parameters  include  radar  characteristics,  radar  placement,  target  geometry/ 
orientation  and  electromagnetic  environmental  considerations  for  example  shadow¬ 
ing  optic.  Other  options  include  range  bins  and  target  compensation  options. 
An  example  of  a  SarTooEm.dat  file  used  during  this  thesis  is  shown  in  Table(A.3.1) 
overleaf. 

A. 3.2  Reflector  Files.  The  reflector  files  are  the  means  whereby  the  operator 
defines  composition  and  placement  of  the  radar  target.  As  mentioned  in  Chapter 
3  the  reflector  files  are  normally  generated  by  the  "EUCLID"  software  and  as  such 
are  not.  readily  amenable  to  manual  construction.  An  example  of  the  inputs  to  a 
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lable  A.l.  SarTool<m  Sensor  and  System  Data  (Sartool.dat) 


exp  1 ' 

TITLE  OF  SIMULATION  RUN 

dihed' 

SPECIFIES  INPUT  REFLECTOR  FILE  I 

'O' 

X  POSITION  OF  SENSOR 

'O' 

V  POSITION  OF  SENSOR 

Z  POSITION  OF  SENSOR 

'00' 

FLIGHT  PATH  ANGLE  OF  SENSOR  ■ 

'O' 

DIVE  ANGLE  OF  SENSOR  j 

'O' 

AVERAGE  SENSOR  SPEED  (SAR  MODE) 

'O' 

X  POSITION  OF  TARGET  REL  TO  MAP  CENTER  j 

'O' 

V  POSITION  OF  TARGET  REL  TO  MAP  CENTER  ! 

‘O' 

Z  POSITION  OF  TARGET  REL  TO  MAP  CENTER 

•u- 

EULER  ANGLE  PHI  FOR  TGT  FRAME  ROTN 

'O' 

EULER  ANGLE  THETA  FOR  TGT  FRAME  ROTN  j 

'O' 

EULER  ANGLE  PSI  FOR  TGT  FRAME  ROTN  1 

TT  GHZ' 

TRANSMITTED  FREQUENCY  ; 

'2.0' 

ELEVATN  3db  BEAMWIDTH  (DEG) 

'2.0' 

AZIMUTH  3db  BEAMWIDTH  (DEG) 

'.or 

ELEVATION  BORESIGHT  ANGLE  (DEG) 

'00.0' 

AZIMUTH  BORESIGHT  ANGLE  (DEG) 

'2.0' 

RANGE  RESOLUTION  (m) 

'2.0' 

AZIMUTH  RESOLUTION  (m) 

T 

NUMBER  OF  RANGE  BINS 

T 

NUMBER  OF  AZIMUTH  BINS 

'200' 

NUMBER  OF  PULSES  TRANSMITTED  j 

T 

TARGET  DATABASE  RECALCULATION  I/VAL 

T 

TERRAIN  GENERATION  OPTION  (SAR) 

'O' 

AVERAGE  TERRAIN  RCS  (dbsm) 

m 

SAR  SPOTLIGHT  OPTION 

V 

ISAR  MODE  OPERATION 

'0' 

ASPECT  ANGLE  INCREMENT  FOR  QKTOOL 

T  MOTION  COMPENSATION  OPTION 

T  REFLECTOR  THINNING  OPTION 


!0’ 

THINNING  THRESHOLD  VALUE  (dbsm) 

Y 

TARGET  SHADOWING  OPTION 

Y 

TERRAIN  SHADOWING  OPTION 

T 

REFLECTOR  I/?  REPORT  OPTION 

H 

REFLECTOR  D/BASE  REPORT  OPTION 

Y 

RANGE  PROFILE  OUTPUT  OPTION 

Table  A. 2.  Reflector  File 


' DM ' 

DIHEDRAL  IDENTIFIER  |! 

'x' 

X  POSITION  OF  DIHED  CENTER  IN  TGT  FRAME 

V 

Y  POSITION  OF  DIHED  CENTER  IN  TGT  FRAME  !j 

V 

Z  POSITION  OF  DIHED  CENTER  IN  TGT  FRAME 

phi’ 

TRANSFORMATION  ANGLE  FROM  REFL  TO  TGT  FRAME  i| 

'theta’ 

TRANSFORMATION  ANGLE  FROM  REFL  TO  TGT  FRAME  :: 

'psi' 

TRANSFORMATION  ANGLE  FROM  REFL  TO  TGT  FRAME  jj 

•dr 

LENGTH  OF  DIHED  IN  CREASE  DIMN  OF  REF  FR.a.mE  H 

I  ’d2’ 

LENGTH  OF  PLATE  PERP  TO  CREASE  DIMN 

{  'sop'' 

SOLID  OBJECTS  PRIMITIVES  INDICES 

reflector  data  file  is  shown  in  Table  A. 2. 

.4. 3. 3  Image. Opts  File.  This  file  specifies  several  processing  and  imaging  op¬ 
tions  used  in  SarTool(m's  imaging  algorithms.  These  options  include  weighting  func¬ 
tions.  scaling  functions  and  the  determination  of  dynamic  ranges.  Because  the  tra¬ 
ditional  imaging  representation  of  SarTooPm  data  was  not  used  in  this  thesis  the 
image. opts  file  will  not  be  discussed  further  although  the  interested  reader  may  find 
additional  information  in  Reference! [26]. 

A. 3. 4  SnrTooltm  Output  Files.  There  are  two  main  output  files  associated 
with  the  above  operator  input  files.  These  two  files  are  briefly  described  below: 

1.  Range  Profile  Data  Base  (RPDB)  file.  This  file  contains  in  phase  and  quadra¬ 
ture  phase  data  for  the  four  polarimetric  components.  This  file  is  created  every 


time  SarTool(m  is  invoked  with  the  associated  data  being  stripped  by  the  db 


strip  program  whl'T  »s  contained  in  Appendix(R). 

2.  Statistics  File.  This  file  contains  the  history  of  the  processed  output.  The 
statistics  file  did  not  play  a  prominent  pait  in  this  thesis  effort  and  will  then-fore 
not  be  discussed  further  in  this  Appendix. 

A. 3.5  Sensor  Geometry  Parameters.  All  of  the  input  files  necessitate  some 
degree  of  translation  between  the  ground  and  target  frames;  accordingly  Figure!  A.l  ) 
demonstrates  the  relationship  between  these  two  co-ordinate  frames. 
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Figure  A. 2.  SarTool(m  Euler  Angles  of  Rotation 


The  Euler  angles  referred  to  in  the  input  files  determine  the  degree  of  rota¬ 
tion  of  the  detailed  axis.  The  actual  rotations  involved  are  described  below  and 
demonstrated  in  Figure  A. 3. 5. 

1.  Euler  Angle  $.  The  angle  $  is  used  to  define  the  amount  by  which  the  x  axis 
is  rotated  towards  the  y  axis  about  the  z  axis. 

2.  Euler  Angle  0.  The  angle  0  is  used  to  define  the  amount  by  which  the  y  axis 
is  rotated  towards  the  z  axis  about  the  x  axis. 

3.  Euler  Angle  The  angle  'F  is  used  to  define  the  amount  by  which  the  x  axis 
is  rotated  towards  the  y  axis  about  the  z  axis. 
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Appendix  B.  Source  Code  of  Thesis  Software 


As  mentioned  in  Chapter  3  this  thesis  eTort  entailed  the  design  of  several 
software  packages  so  tha’  polarimetric  behaviour  could  be  readily  examined.  The 
operational  function  of  this  software  is  described  in  Chapter  3.  To  assist  in  any  future 
research  endeavours  in  this  field  a  source  listing  of  each  of  the  packages  appears  below: 

B.l  Reflector  File  Software 

The  following  piece  of  software  belongs  to  the  sponsor  of  this  thesis  effort  and 
was  written  by  M.  Bryant. 


/***** 

M0D_RFL.C  -  Program  to  enable  easy  modification  to  SarTool 
formatted  reflector  database  files. 

*****/ 


#include  <stdio.h> 

#mclude  <string.h> 

^include  <math.h> 

ttdefine  max_primitives  100 

PILE  *dat_file,  *input_file,  *output_f ile ; 
char  new.name [80] ,  temp [80]  ,  header [80] ; 

int  num.FP,  num_DH,  num.TH,  num.EL,  num_CD,  num_CY,  num_TP, 
num.EP,  num_CC,  num_CV,  num_UG,  num_CN,  num_ET,  index; 

typedef  struct 
{  double  x,  y,  z; 

}  cartesian; 
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cartesian  start,  end,  increment; 


typedef  struct 
{  double  phi,  theta,  psi; 

}  euler; 

typedef  struct 
{  char  type [3] ; 
cartesian  location; 
euler  orientation; 
union 

{  struct  {  double  xl,  x2,  y2,  x3,  y3;  >  FP; 
struct  {  double  dl,  d2 ;  }  DH; 
struct  {  double  d,  mxpkdb ;  }  TH; 
struct  {  double  axisa,  axisb,  axisc;  }  EL; 
struct  {  double  dl,  d2,  rcurv,  pkdbm;  >  CD; 
struct  {  double  axisa,  axisb,  height;  >  CY ; 
struct  {  double  axisa,  axisb,  height;  >  TP; 
struct  {  double  a,  b;  }  EP; 

struct  {  double  a,  b,  height,  xbmvth,  ybmvth,  pkdbsm;  > 


CC; 

struct  {  double  a,  b,  height,  xbmwth,  ybmwth,  pkdbsm;  > 

Ob¬ 
struct  <  double  a,  length;  >  OG; 
struct  {  double  a,  b,  lenl,  len2;  }  CN; 
struct  {  double  axisa,  axisb,  cyoffx, 

angrot,  axisc,  axisd,  height;  }  £T ; 

}  unique; 

int  sopa,  sopb,  sopc; 

}  primitive; 

primitive  ptmax.primitives] ; 
int  i ; 

void  read_header() 

{  if  (NULL**fgets (header,  80,  input.f lie)) 

{  puts("*****Error  reading  header*****\n") ; 
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exit  ()  ; 

} 

} 


void  read.FP (num) 
int  num; 

{  ii  ( (13 ! =fscanf (input .file , 


"*/,s  ‘/.If  ’/.If  ’/.If  */.lf  ‘/.If  ‘/.If  ‘/.If 


‘/.If  '/.If  ‘/.If  '/.If  y.d\n", 
pCi] .type, 

&p [ i] .location. x,  &p [1] . location . y ,  ftp [i]  . location. z , 


&p[i] .orientation. phi,  &p[i] .orientation. theta, 

&p[i] . orientation. psi, 

&p [i] . unique. FP.xl,  &p[i] . unique. FP.x2, 

&p[i] . unique. FP.y2, 

&p[i] . unique. FP.x3,  &p[i] . unique. FP.y3, 

&p [i] .sopa)) I  I (0!=strcmp(p[i] .type,  "FP"))) 

{  printf ("*****Error  reading  flatplate  '/,d*****\n" ,  num) ; 


exit()  ; 

> 

++i; 

> 

void  read_f latplatesO 
{  int  num; 

if  (1 ! =f scanf (input_f ile,  "'/,d\n",  ftnum.FP)) 

{  puts("*****Error  reading  number  of  f latplates*****\n") ; 


exit()  ; 

> 

for  (num=0 ;num<num_FP;++num)  read_FP(nnm) ; 

> 

void  read_DH(num) 
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int  num; 

{  if  ((11 !=fscanf (input.file,  "s/,s  '/.If  ‘/.If  */,lf  */,lf  '/.If  '/.If  '/.If 


•/.If  y.d  */,d\n"  , 
p[i]  .type, 

&p[i] .location. x,  &p [ i] .location. y,  ftp [i] .location. z , 


&p[i] . orientation .phi ,  &p[i] .orientation. theta, 

&p[ij . orientation. psi, 

ftpCi] . unique. DH.dl,  &p[i] . unique . DH . d2 , 

&p[i].sopa,  ftp [1] . sopb) ) | | (0 ! =strcmp(p [i] .type ,  ”DH”))) 


{  printf ("*****Error  reading  dihedral  ’/,d*****\n",  num); 


exit  ()  ; 

} 

++i; 

} 

void  read_dihedrals() 

{  int  num; 

if  (1 ! =f  scanf  (input_f  ile,  "'/.d\n"  ,  ftnum.DH)) 

{  puts("*****Error  reading  number  of  dihedrals*****\n") ; 


exitQ  ; 

> 

for  (num  =  0 ;num<num_DH;++num)  read.DH(num) ; 

} 

void  read_TH(num) 
int  num; 

{  if  (( 12  !=f  scanf  (input. file,  "'/.s  %lf  ’/.If  ’/.If  ‘/.If  ’/.If  ’/.If  ’/.If 

•/.if  */.d  y.d  y.d\n". 
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P  [i] -type, 

&p[i] .location. x,  ftp[i] .location. y,  ftp[i] .location. z, 


&p[i] .orientation. phi,  &p[i] . orientation. theta, 

&p [i] . orientation .psi , 

&p[i] . unique. TH.d,  &p[i] .unique .TH.mxpkdb, 

&p[i].sopa,  &p[i].sopb, 

&p[i] .sopc)) I  I (0 !=strcmp(p[i] .type,  "TH"))) 

{  printf ("*****Error  reading  trihedral  '/,d*****\n",  num); 


exit()  ; 

> 

++i; 

} 

void  read_trihedrals() 

{  int  num; 

if  (1 !=fscanf (input_file,  ’’’/.dXn",  ftnum_TH)) 

{  puts("*****Error  reading  number  of  trihedrals*****\n”) ; 


exit() ; 

} 

for  (num  =  0 ;nura<num_TH;++num)  read_TH(num) ; 

} 

void  read_EL(num) 
int  num; 

{  if  ((11  !=fscaaf (input_file,  "%s  '/,lf  '/.If  '/.If  '/.If  Ilf  '/.If  '/.If 


’/.if  y.if  y.d\n", 

p[i] -type, 

ftp [i] . location .x ,  ftp[i] .location. y,  ftp[i] .location. z, 


ftp[i] .orientation. phi,  ftp [i] .orientation. theta, 
&p[i] .orientation. psi, 
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&p[i] . unique . EL . axisa,  &p[i] . unique . EL . axi sb , 

&p[i] .unique. EL. axisc,  &p[i] .sopa)) 

I  I (0!=strcmp(p[i] .type,  "EL"))) 

{  printf  ("*****Error  reading  Ellipsoid  '/,d*****\n"  ,  num)  ; 


exit()  ; 

> 

++i; 

} 

void  read_ellipsoids () 

{  mt  num; 

if  (1  !=fscanf (input_f ile ,  "'/,d\n",  &num_EL)) 

{  puts("*****Error  reading  number  of  ellipsoids*****\n") ; 


exit  ()  ; 

} 

for  (num  =  0 ; num<num_EL ; ++num)  read.EL(num) ; 

} 

void  read_CD(num) 
int  num; 

{  if  (( 13 !=fscanf (input .file,  "'/.s  Xlf  Xlf  '/.If  7.1f  ‘/.If  ’/.If  ’/.If 


'/.If  '/.If  Xlf  Xd  y,d\n" , 
p[i] .type, 

&p[i] .location. x,  &p[i] .location. y,  4p[i] .location. z, 


&p[i] .orientation. phi,  ftp[i] .orientation. theta, 

&p[i] . orientation. psi, 

&p[i] . unique. CD. dl ,  4p[i] .unique. CD. d2, 

&p[i] .unique. CD. rcurv,  4p[i] .unique .CD .pkdbm, 

&p[i].sopa,  &p[i].sopb)) 

I  I (0 ! =strcmp(p[i] .type,  "CD"))) 

{  printf  ("*****Error  reading  Curved  Dihedral  Xd*****\n" , 
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num)  , 
exit ()  ; 

> 

++i; 

> 

void  read_curveddihedrals() 

{  int  num; 

if  (1  !=fscanf (input _file,  "*/,d\n",  &num_CD)) 

{  puts("*****Error  reading  number  of  curved 
dihedrals*****\n") ; 
exit  ()  ; 

> 

for  (num  =  0 ;num<num_CD ;++num)  read_CD(num) ; 

> 

void  read_CY(num) 
int  num; 

{  if  ((11!  =f scanf  ( input _file ,  "*/.s  ‘/.If  ‘/.If  */.lf  ‘/.If  */.lf  */.lf  */.lf 

*/.lf  */.lf  '/.dXn"  , 

P  [i] -type, 

&p [i] . location .x ,  &p[i] .location. y,  &p[i] .location. z, 


&p[i] .orientation. phi,  &p[i] .orientation. theta, 

&p[i] . orientation. psi, 

&p Ci] . unique. CY.axisa,  &p[i] . unique. CY.axisb, 

&p[i] . unique. CY. height,  &p[i] .sopa)) 

I  I (0!=strcmp(p[i] .type,  "CY"))) 

{  printf ("*****Error  reading  Cylinder  %d*****\n",  num); 


exitO  ; 

> 

++i ; 

> 

void  read_cylinders() 
{  int  num; 
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if  (1  !=fscanf (input_f ile ,  "'/.d\n",  &num_CY)) 

{  puts ("*****Error  reading  number  of  cylinders*****\n") ; 


exit  ()  ; 

> 

for  (num  =  0;num<num_CY;++num)  read_CY (num) ; 

> 

void  read_TP(num) 
int  num; 

{  if  ( (12  !  =fscanf  (input _f  ile,  '7.s  '/.If  ‘/.If  ‘/.If  ‘/.If  ‘/.If  '/.If  '/.If 

‘/.If  ‘/.If  ‘/.d  */,d\n"  , 
p[i] .type, 

&p [i] . location. x,  &p[i] .location. y,  &p[i] .location. z. 


ftp[i] .orientation. phi,  &p[i] .orientation. theta, 

&p  Ci3 . orientation. psi, 

&p[i] . unique. TP. axisa,  &p[i] .unique. TP. axisb, 

&p[i] .unique .TP .height ,  &p[i] .sopa, 

&p[i] . sopb) ) 

I  I (0!=strcmp(p[i] .type,  "TP"))) 

{  printf  ("*****Error  reading  Top  Hat  y,d*****\n" ,  num); 

exit  () ; 

> 

++i; 

> 

void  read_tophats() 

{  int  num; 

if  (1 ! =f scanf (input.f ile,  "'/,d\n"  ,  &num_TP)) 

{  puts("*****Error  reading  number  of  top  hats*****\n") ; 


exit() ; 

> 

for  (num  *  0 ; num<num_TP ; ++num)  read_TP(num) ; 
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} 


void  read_EP(num) 
in*  num; 

{  if  (( 10  !=fscanf  (input  .file,  ’7.s  '/.If  ‘/.If  ‘/.If  */.lf  */.lf  ‘/.If  '/.If 

’/.If  */.d\n"  , 
p  [1] . type , 

&p [1] . location .x ,  &p[i] . location. y,  &p  [i] . location . z , 


&p[i] . orientat ion . phi ,  fep[i] .orientation. theta, 

&p[i] . orientation .psi , 

&p[i] . unique. EP. a,  &p[i] . unique. EP.b, 

&p [ij . sopa) ) 

11(0' =strcmp(p [1] . type ,  "EP") ) ) 

{  prmtf  ("*****Error  reading  Elliptical  Plate  '/,d*****\n" , 


num)  ; 
exitO  ; 

} 

++i; 

} 

void  read_ellipticalplates() 

{  int  num; 

if  (1  !=fscanf (input_file,  "'/,d\n",  &num_EP)) 

{  puts("*****Error  reading  number  of  elliptical 
plates*****\n") ; 
exitO  ; 

> 

for  (num  *  0 ;num<num_EP;++num)  read_EP(num) ; 

> 

void  read_CC(num) 
int  num; 

{  if  (( 14  !«fscanf  (input,  file,  "'/.s  %lf  ‘/.If  Vf  Xlf  '/.If  ’/.If  '/.If 

•/.if  ‘/.if  y.if  */.if  y.if  y.d\n" , 
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p[i] -type, 

&p[i] .location. x,  ftp  [  i]  . location . y ,  ftp  [1] . location . z , 


&p[i] . orientat ion . phi ,  ftp[i] . orientation . theta , 
ftp  C i ] . orientation .psi , 

&p[i] . unique. CC. a,  ftp [i] . unique .CC . b , 

&p[i] .unique.CC. height,  ftp [i] .unique .CC.xbmwth, 
ftp[i] . unique. CC.ybmwth,  ftp[i] . unique . CC . pkdbsm , 
ftp [1] . sopa) ) 

I | (0 ! =strcmp(p [i] .type,  "CC"))) 

{  printf ("**** *Error  reading  Complex  Cavity  '/,d*****\n"  , 

num)  ; 
exit  ()  ; 

> 

++i; 

> 

void  read_complexcavities() 

{  int  num; 

if  (1 !  =f scanf (input.f ile,  "5£d\n",  ftnum_.CC)) 

{  puts ("*****Error  reading  number  of  complex 
cavities*****\n") ; 
exit ()  ; 

} 

for  (num  =  0 ; num<nura_CC ; ++num)  read_CC(num) ; 

> 

void  read_CV(num) 
int  num; 

{  if  (( 14 !  =fscanf  (input _file ,  "*/,s  y,lf  ‘/.If  '/.If  '/.If  '/.If  */«lf 


'/.If  ‘/.If  ‘/.If  '/.If  '/.If  '/. d\n", 
p [i] -type, 

ftp[i] .location. x,  &p[i] .location. y,  ftp [i] .location. z, 


ftp Ci] .orientation. phi,  ftp Ci] . orientation. theta , 

&p[i] . orientation .psi , 

ftp [i] . unique. CV. a,  ftp  Ci] .unique .CV.b, 
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&p[i] . unique. CV. height ,  &p[i] . unique. CV.xbmwth, 
ftp[i] .unique .CV.ybmwth,  ftp[i] .unique .CV .pkdbsm , 
ftp  [i]  .  sopa) ) 

I  I (0 ! =strcmp(p[i] .type,  "CV"))) 

{  prmtf  ("*****Error  reading  Cavity  '/,d*****\n"  ,  num)  ; 
exit ()  ; 

> 

++i; 

> 

void  read_cavities() 

{  int  num; 

if  ( 1 !  =f scanf (input_f lie ,  "*/,d\n",  ftnum_CV)) 

{  puts("*****Error  reading  number  of  cavities*****\n") ; 


exit  ()  ; 

> 

for  (num  =  0 ; num<num_CV ; ++num)  read_CV(num) ; 

> 

void  read_0G(num) 
int  num; 

{  if  (( 10  !-f  scanf  (input  .file,  "*/.s  ’/.If  ’/.If  '/.If  '/.If  '/.If  '/.If  */.lf 

'/.If  Xd\n", 

p[i] .type, 

&p[i] .location. x,  ftp  [i] . location. y,  ftp [i] . location .z , 


ftpLi] .orientation. phi,  &p[i] .orientation. theta, 

&p[i] . orientation. psi, 

&p[i] . unique. OG. a,  &p[i] . unique . 0G . length , 

&p[i] .sopa)) 

I  I (0!=strcmp(p[i] .type,  "OG"))) 

{  printf ("*****Error  reading  Ogive  '/,d*****\n" ,  num); 
exit ()  ; 

> 

++i; 

> 
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void  read_ogives() 

{  mt  num; 

if  ( 1 !  =f scanf (input_f ile ,  "’/dXn",  &num_0G)) 

{  puts("*****Error  reading  number  of  ogives*****\n") ; 
exit ()  ; 

> 

for  (num  =  0 ; num<num_QG ; ++num)  read_0G (num) ; 

} 

void  read_CN(num) 
int  num; 

{  if  ( (12 !  =f  scanf  (input  _f  ile ,  "7.s  */,lf  ’/.If  ’/.If  ’/.If  ’/.If  ’/.If  ’/.If 

'/.If  ’/.If  ’/.If  */.d\n" , 
p  [i] .type, 

&p [ i] .location. x,  &p[i] .location. y,  &p[i] .location. z, 


&p[i] .orientation. phi,  &p[i] . orientation. theta, 

&p[i] . orientation .psi , 

&p[i] . unique. CN. a,  &p[i] . unique. CN.b, 

&p[i] . unique. CN.lenl,  &p[i] . unique. CN.len2, 

&p[i] .sopa)) 

11(0! -strcmp(p[i] .type,  "CN") ) ) 

{  printf  ("****  *Error  reading  Cone  '/,d*****\n"  ,  num); 
exit () ; 

> 

++i; 

> 

void  read_cones() 

{  int  num; 

if  (1  !=fscanf (input_file,  "*/,d\n",  4num_CN)) 

{  puts(" :'****Error  reading  number  of  cones*****\n") ; 
exitQ  ; 

> 

for  (num  -  0;num<num_CN;++num)  read_CN(num) ; 

> 

void  read_ET(nu m) 
int  num; 
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{  if  ( (17  '=fscanf  (input  .file,  "V.s  */.lf  ’/.If  */.lf  ‘/.If  %lf  ’/.If  ’/.If 
*/.if  ’/.If  ‘/.if  ’/.If  ’/.If  V.lf  */.lf  */.d  '/.d\n" , 

p  Ci] -type , 

&p[i] .location. x,  &p [1] . location. y,  &p[i] .location. z, 


&p[i] .orientation. phi,  &p[i] .orientation. theta, 

&p[i] . orientation. psi, 

&p[i] .unique . ET. axisa,  &p[i] . unique . ET . ax isb , 

&p[i] . unique. ET.cyoffx,  &p[i] . unique. ET.cyoffy, 

&p[i] . unique. ET.angrot,  &p[i] . unique. ET.axisc, 

&p[i] . unique. ET.axisd,  &p[i] . unique . ET. height , 
&p[i].sopa,  &p[i].sopb)) 

I  I (0 !=strcmp(p[i] .type,  "ET"))) 

{  printf ("*****Error  reading  Elliptical  Base  Top  Hat 
'/,d*****\n" ,  num)  ; 
exit() ; 

> 

++i; 

> 

void  read_ellipticalbasetophats () 

{  int  num; 

if  (1  !=fscanf (input.file,  "'/,d\n",  &nuo_ET)) 

{  puts("*****Error  reading  number  of  elliptical  base  top 


hats*****\n") ; 
exit()  ; 

> 

for  (num  =  0;num<num_ET;++num)  read_ET(num) ; 

> 

void  write_header() 

{  if  (E0F==fprintf (output.f ile,  "%s\n",  header)) 
{  puts ("*****Error  writing  header*****\n") ; 
exit()  ; 

> 

> 
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void  write_FP(num) 
int  num; 

{  if  (EOF==fprintf (output_f ile,  “5! 'As  '/,10.31f  */,10.31f  '/10.31f 

*/,10 . 31f  */,10.31f  7,10 . 31f  '/,10 . 31f  */.10.31f  */.10.31f  */.10.31f  Xl0.31f 

*/,5d\n"  , 
p[i] -type, 

p[i] .location. x,  p[i] .location. y,  p[i] . location. z, 
p[i] . orientation .phi ,  p[i] .orientation. theta, 
p[i] . orientation . psi , 

p [1] .unique . FP .xl ,  p [i] . unique . FP . x2 ,  p [i] .unique . FP ,y2 , 


p[i] . unique. FP.x3,  p[i] . unique. FP .y3, 
p [i] .sopa)) 

{  printf ("*****Error  writing  flatplate  '/,d*****\n" ,  num); 


exit  ()  ; 

> 

++i; 

} 

void  write.f latplatesO 
{  int  num; 

if  (E0F==f printf (output _f ile ,  "%5d\n" ,  num.FP)) 

{  puts("*****Error  writing  number  of  flatplates*****\n") ; 


exit() ; 

> 

for  (num  *  0;num<num_FP;++num) 
{  write_FP(num) ; 

} 

> 

void  write_DH(num) 
int  num; 
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{  if  (EQF==fprintf (output.f lie ,  "*!As  '/,10.31f  '/,10.31f  '/,10.31f 

'/,10.31f  */,10.31f  */,10 . 31f  y,10.31f  */,i0 . 31f  */,5d  '/.5d\n", 
p [i] .type, 

p[i] .location. x,  p[i] .location. y,  p [i] . location .z , 
p[i] .orientation. phi,  p[i] .orientation. theta, 
p[i] . orientation . psi , 
p[i] . unique. DH.dl ,  p[i] . unique. DH.d2, 
p[i].sopa,  p  [i]  .  sopb) ) 

{  printf  ("*****Error  writing  dihedral  ‘/,d*****\n" ,  num) ; 

exit()  ; 

> 

++i; 

> 

void  write_dihedrals() 

{  int  num; 

if  (E0F==f printf  (output.file ,  "'/,5d\n"  ,  num.DH)) 

{  puts("*****Error  writing  number  of  dihedrals*****\n") ; 


exit() ; 

} 

for  (num  =  0 ;num<num_DH;++num)  write_DH(num) ; 

> 

void  write_TH(num) 
int  num; 

{  if  (EOF==fprintf (outpuc.f ile,  "'/,4s  '/,10.31f  '/,10.31f  ’/.10.31f 

'/.10 . 31f  y,10.31f  y,10.31f  '/,10.31f  y,10.31f  y,5d  y,5d  y,5d\n", 
p[i] -type, 

p[i] .location. x,  p[i] .location. y,  p[i] .location. z, 
p[i] .orientation. phi,  p[i] .orientation. theta, 
p[i] .orientation. psi, 

p[i] .unique .TH.d,  p[i] . unique. TH.mxpkdb, 
p[i] .sopa,  p[i] .sopb,  p[i].sopc)) 

{  printf ("*****Error  writing  flatplate  y,d*****\n",  num); 
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exit  ()  ; 

} 

++i; 

> 

void  write.tnhedralsQ 
{  int  num; 

if  (EOF==fprintf (output .file,  "*/,5d\n",  num.TH)) 

{  puts("*****Error  writing  number  of  f latplates*****\n") ; 


exit  () ; 

> 

for  (num  =  0;num<num_lH;++num)  write_FP(num) ; 


void  write_EL(num) 
int  num; 

{  if  (EOF==fprintf (output .file,  "'/As  '/,10.31f  %10.31f  Xl0.31f 

*/,10 . 31f  y.l0.31f  */.10.31f  */,10 . 31f  7.10. 31f  */.10.31f  */.5d\n'', 
pti] -type, 

p[i] .location. x,  p[i] .location. y.  p[i] .location. z, 
p C i] . orientation. phi ,  p[i] .orientation. theta, 
p[i] . orientation. psi, 

p [ i] .unique. EL. axisa,  p[i] .unique. EL. axisb, 
p[i] .unique. EL. axisc,  p[i].sopa)) 

{  printf ("*****Error  writing  Ellipsoid  y,d*****\n",  num); 


exit  () ; 

> 

++i; 

> 

void  write_ellipsoids() 
{  int  num; 
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if  (EOF==f printr  (output.f ile ,  "75d\n",  num_EL)) 

{  puts("*****Error  writing  number  of  ellipsoids*****\n") ; 


exit()  ; 

> 

for  (num  =  0;num<num_EL;++num)  write_EL(num) ; 

> 

void  write_CD(num) 
int  num; 

{  if  (EOF==fprintf (output _f ile,  "74s  7,10. 31f  7,10.31f  '/,10.31f 

'/.10 . 31f  7,10 . 31f  7,10. 31f  7.10. 31f  7,10. 31f  7,10. 31f  7,10. 31f  7.5d  7.5d\n" 


p[i] -type, 

p[i] .location. x,  p[i] .location. y,  p[i] .location. z, 
p[i] .orientation. phi,  p[i] .orientation. theta, 
p[i] .orientation.psi, 
p[i] . unique. CD. dl,  p[i] .unique. CD. d2, 
p[i] . unique. CD .rcurv,  p[i] . unique. CD. pkdbm, 
p[i].sopa,  p  [i]  .  sopb)  ) 

{  printf  ("*****Error  writing  Curved  Dihedral  7.d*****\n" , 


num) ; 
exit  ()  ; 

> 

++i; 

> 

void  write_curveddihedrals() 

{  int  num; 

if  (E0F==fprintf (output.f ile,  "7,5d\n",  num.CD)) 
{  puts("*****Error  writing  number  of  curved 
dihedrals*****\n") ; 
exitO  ; 

> 

for  (num  *  0 ;num<num_CD;++num)  write_CD(num) ; 


} 

void  write_CY(num) 
int  num; 

{  if  (EOF==fprintf (output.f ile,  "'/,4s  Xl0.31f  %10.31f  '/I0.31f 

*/,10 . 31f  */.10.31f  */.10.31f  '/,10 . 31f  y,10.31f  */.10.31f  '/.5d\n"  , 
p  [i] -type, 

p[i] .location. x,  p[i] .location. y,  p [i] .location. z, 
p[i] .orientation. phi,  p[i] .orientation. theta, 
p[i] . orientation. psi , 

p[i] . unique. CY.axisa,  p[i] . unique. CY.axisb, 
p[i] .unique .CY .height ,  pti] .sopa)) 

{  printf ("*****Errcr  writing  Cylinder  '/,d*****\n" ,  num); 

exit ()  ; 

> 

++i; 

> 

void  write_cylinders() 

{  int  num; 

if  (E0F==f printf  (output  .file,  '"/,5d\n"  ,  num_CY)) 

{  puts("*****Error  writing  number  of  cylinders*****\n") ; 


exitO  ; 

> 

for  (num  =  0 ;num<num_CY;++num)  write_CY(num) ; 

> 

void  write_TP(num) 
int  num; 

{  if  (EOF«fprintf (output_f ile,  "'/,4s  y,10.31f  y,10.31f  '/,10.31f 

y,10.31f  y,10.31f  '/,10.31f  y,10.31f  y,10.31f  y,10.31f  y,5d  y.5d\n", 
p  [i] .type, 

p[i] . location. x,  p[i] .location. y,  p[i] .location. z, 
p[i] .orientation. phi,  p[i] .orientation. theta, 
p[i] . orientation. psi , 
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p[i] .unique. TP. axisa,  p[i] . unique. TP. axisb, 
p[i] . unique . TP .height ,  p[i] .sopa, 
p [i] . sopb) ) 

{  printf  ("*****Error  writing  Top  Hat  Xd*****\n" ,  num); 

exit  ()  ; 

> 

++i; 

> 

void  write_tophats() 

{  int  num; 

if  (EOF==fprintf  (output _file  ,  "'/.5d\n"  ,  num_TP)) 

{  puts("*****Error  writing  number  of  top  hats*****\n") ; 

exit ()  ; 

> 

for  (num  =  0;num<num_TP;++num)  write  TP(num); 

> 

void  write.EP(num) 
int  num; 

{  if  (EOF==fprintf (output .file,  "74s  Xl0.31f  y.l0.31f  '/.10.31f 

«10 . 31f  7.10. 31f  */,10.31f  */.10.31f  %10 . 31f  ’/,5d\n"  , 
p[i]  -type, 

p[i] .location. x,  p[i] .location. y,  p[i] .location. z, 

p[i] .orientation. phi ,  p[i] . orientation .theta, 

p  C i] . orientation. psi, 

p[i] . unique. EP. a,  p[i] . unique. EP.b, 

pCi] .sopa)) 

{  printf  ("*****Error  writing  Elliptical  Plate  '/,d*****\n" , 

num) ; 
exit()  ; 

> 

♦+i; 

> 

void  write.ellipticalplatesQ 
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{  int  num; 

if  (EQF==fprintf (output_f ile ,  "'/,5d\n",  num_EP)) 

{  puts("*****Error  writing  number  of  elliptical 
plates*****\n") ; 
exit() ; 

> 

for  (num  =  0;num<num_EP;++num)  write_EP(num) ; 

> 

void  write_CC(num) 
int  num; 

{  if  (EOF==fprintf (output.f ile,  "%4s  */,10.31f  */,10.31f  '/, i0.31f 

’/.10.31f  */,10 . 31f  '/,10 . 31f  y,10.31f  */,10 . 31f  */.10.31f  */.10.31f  */.10.31f 

*/,10 . 31f  '/,5d\n"  , 
p[i] .type, 

&p[i] .location. x,  &p [i] . location . y ,  &p[i] .location. z, 

&p[i] .orientation. phi,  &p[i] . orientation. theta, 

&p[i] . orientation. psi, 

&p[i] . unique. CC. a,  &p[i] . unique. CC.b, 

&p[i] . unique. CC. height ,  &p[i] .unique .CC .xbmwth, 

&p[i] .unique .CC.ybmwth,  &p[i] . unique. CC.pkdbsm, 

&p[i] .sopa)) 

{  printf  ("♦♦♦♦♦Error  reading  Complex  Cavity  '/,d*****\n" , 

num) ; 
exit ()  ; 

> 

++i; 

} 

void  write_complexcavities() 

{  int  num; 

if  (EOF=*fprintf (output.f ile ,  "'/,5d\n”,  num_CC)) 

{  puts("*****Error  writing  number  of  complex 
cavities*****\n") ; 
exit()  ; 

> 

for  (num  *  0 ; num<num_CC ; ++nura)  write_CC(num) ; 
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> 

void  write_CV(num) 
int  num; 

{  if  (EOF==fpnntf (output _f lie,  "'/,4s  '/,10.31f  '/,10.31f  '/,10.31f 

'/.10.31f  '/,10 . 31f  */.10.31f  '/,10 . 31f  */,10 . 31f  */.10.31f  '/.10.31f  */.10.31f 

'/.10.31f  '/.5d\n”, 
p[i] -type, 

p[i] .location. x,  p[i] .location. y,  p[i] .location. z, 

pCi] .orientation. phi,  p[i] .orientation. theta, 

p[i] . orientation. psi, 

p[i] . unique. CV. a,  p[i] . unique. CV.b, 

p[i] . unique. CV. height,  p[i] . unique. CV.xbmwth, 

p[i] . unique. CV.ybmwth,  p[i] . unique. CV.pkdbsm, 

p[i] .sopa)) 

{  printf ("*****Error  writing  Cavity  y,d*****\n",  num); 
exit()  ; 

> 

++i; 

} 

void  write_cavities() 

{  int  num; 

if  (EQF==fprintf (output .file,  "'/,5d\n"  ,  num.CV)) 

{  puts("*****Error  writing  number  of  cavities*****\n") ; 

exitQ  ; 

} 

for  (num  =  0 ;num<num_CV;++num)  write.CV(num) ; 

> 

void  write_OG(num) 
int  num; 

{  if  (E0F==f  printf  (output  .file,  "y.4s  y,10.31f  y,10.31f  '/.lO.Slf 

'/.10 . 31f  y,10.31f  y,10.31f  y,10.31f  y,10.31f  y,5d\n”, 

p[i] .type, 

p[i] .location. x,  pCi] .location. y,  p[i] .location. z, 
p[i] .orientation. phi,  pCi] .orientation. theta, 
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p[i] . orientation . psi , 

p[i] . unique . OG . a,  p[i] . unique. OG. length , 
p[i]  .  sopa)) 

{  printf  ("*****Error  writing  Ogive  ‘/,d*****\n" ,  num)  ; 
exit  ()  ; 

} 

++i; 

> 

void  write_ogives () 

{  int  num; 

if  (EQF==fprintf (output_f ile ,  "'/.SdXn1',  num_0G)) 

{  puts ("*****Error  writing  number  of  ogives*****\n") ; 
exit ()  ; 

> 

for  (num  =  0 ;num<num_0G ;++num)  write_OG(num) ; 

} 

void  write_CN(num) 
int  num; 

{  if  (EOF==fprintf ( output _f ile ,  "‘/.4s  '/.10.31f  */.10.31f  ’/.10.31f 
'/,  10 . 31f  y,10.31f  */,10 . 31f  ‘/.10 . 31f  '/,10 . 31f  %10.31f  */.10.31f  */.5d\n'' 


p[i]  .type, 

p[i] .location. x,  p [i] . location . y,  p[i] .location. z, 

p[i] . orientation. phi ,  p[i] . orientation . theta, 

p[i] . orientation .psi , 

p [ i] . unique. CN. a,  p[i] . unique. CN.b, 

p[i] . unique. CN.lenl,  p[i] . unique. CN.len2, 

p[i] .sopa)) 

{  printf  ("*****Error  writing  Cone  */,d*****\n" ,  num)  ; 
exit() ; 

> 

++i; 

} 

void  write_cones() 

{  int  num; 

if  (EOF**fprintf (output.f ile,  ”'/(5d\n",  num.CN)) 

{  puts("*****Error  writing  number  of  cones*****\n") ; 
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exit  ()  ; 

> 

for  (num  =  0 ; num<num_CN ; ++num)  write_CN(num) ; 

> 

void  write_ET(num) 
int  num; 

{  if  (EOF==fprintf (output _f lie ,  "'/.4s  '/,10.31f  '/,10.31f  */,10.31f 

'/.  10 . 31f  '/,10 . 31f  */,10 . 31f  '/,  10 . 31f  '/,10 . 31f  '/.10.31f  '/,10.31f  '/.10.31f 

*/,10 . 31f  '/.10.31f  '/,10 . 31f  ‘/,5d  '/,5d\n"  , 
p[i] .type, 

p  [i] . location. x ,  p  [i] . location . y ,  p [i] . location .2 , 
p[i] .orientation. phi ,  p[i] . or lent at  ion . theta , 
p[i] . orientation . psi , 

p[i] . unique. ET.axisa,  p[i] . unique. ET. axisb, 
p[i] .unique. ET.cyoffx,  p[i]  . unique . ET . cyoffy , 
p[i] . unique . ET . angrot ,  p[i] . unique . ET . axisc , 
p[i] . unique. ET.axisd,  p[i] .unique .ET. height , 
p[i]  .sopa,  p[i]  .  sopb) ) 

{  printf ("*****Error  writing  Elliptical  Base  Top  Hat 
’/.d*****\n" ,  num)  ; 
ex  it  ( )  ; 

} 

++i; 

> 

void  write_ellipticalbasetophats() 

{  int  num; 

if  (E0F==f  prmtf  (output_f  lie  ,  "'/,5d\n"  ,  num_ET)) 

{  puts("*****Error  writing  number  of  elliptical  base  top 


hats*****\n") ; 
exit ()  ; 

> 

for  (num  *  0 ;num<num_ET;++num)  write_ET(num) ; 

> 

void  read_input_f ile() 

{  mput.file  =  fopen("INPUT.RFL",  "r"); 
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if  (input.file  ”=  NULL) 

{ 

puts  ("***  Cfu't  open  input  file  ***"); 
exit  (0) ; 

> 

1=0; 

read_header() ; 
read_f latplatesO  ; 
read_dihedrals () ; 
read_tnhedrals()  ; 
read_ellipsoids() ; 
read.curveddihedralsQ ; 
read_cylmders ()  ; 
read_tophats () ; 
read_ellipticalplates() ; 
read_complexcavities () ; 
read_cavities() ; 
read_ogives () ; 
read_cones () ; 

read_ellipticalbasetophats() ; 

> 

void  vrite_rf l_f ile() 

{  i=0; 

write.header () ; 
vnte_f  latplatesO  ; 
write_dihedrals() ; 
write_tnhedrals() ; 
write_ellipsoids() ; 
wnte_curveddihedrals() ; 
wnte_cylinders() ; 
vnte_tophats()  ; 
wnte.ellipticalplatesO ; 
wnte_complexcavities() ; 
vrite_cavities() ; 
write_ogives() ; 
write_cones() ; 

wnte_ellipticalbasetophats() ; 

> 

cartesian  assign_cartesian(argl ,  arg2 ,  arg3) 
double  argl,  arg2,  arg3; 


B-24 


{ 

cartesian  result; 
result . x=argl ; 
result . y=arg2 ; 
result . z=arg3 ; 
return  (result ) ; 

> 

cartesian  add_cartesian(cl ,  c2) 
cartesian  cl ,  c2  ; 

{ 

cartesian  result; 
result . x=cl . x+c2 . x ; 
result .  y-cl . y+c2 . y ; 
result . z=cl . z+c2 . z ; 
return  (result) ; 

} 

void  translate(index) 
int  index; 

{ 

p [index] . location=add_cartesian(p [index] . location, 


> 

void  calculate_inc() 

{ 

increment . x=(end . x-start . x)/9 . 0 ; 
increment . y= (end . y-start . y) /9 . 0 ; 
increment .  z=  (end . z-s  tai.  L .  z) /9 . 0  ; 

> 

int  modO 

{ 

increment =as sign. cartesian (start) ; 
translateO  ; 
calculate_inc() ; 
for  (i=0 ; i<10 ; ++i) 

{  strcpy(new_name,  "test_fmt"); 
sprintf  (temp ,  "'/,-3d"  ,  i) ; 


increment)  ; 
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strcat (strcat (new.name ,  temp); 

output_file  =  f open(new_name ,  "w"); 
if  (output.file  ==  NULL) 

{  puts  ("***  Can’t  open  output  file  ***"); 
exit  (0)  ; 

> 

write_rf l_f ile() ; 
close(output_f ile) ; 
translateO  ; 

> 

> 

void  read_dat_f ile() 

{ 

dat.file  =  f open("M0D_RFL .DAT",  "r") ; 
if  (input_file  ==  NULL) 

{ 

puts  ("***  Can’t  open  data  file  ***”); 
exit  (0)  ; 

> 

if  (3  !-(fscanf  (dat.file,  "Ilf  '/.If  ‘/.lf\n", 

ftstart.x,  ftstart.y,  istart.z))) 

{  puts("***  Error  reading  data  file  ***"); 
exit (0)  ; 

} 

if  (3 !  =  (fscanf (dat.file,  "'/.If  */,lf  %lf\n",  ftend.x,  &end.y, 
&end.z))) 

{  puts("***  Error  reading  data  file  ***"); 
exit (0)  ; 

> 

if  (1 ! =(f scanf (dat.f ile ,  ”'/,d"  ,  ftindex))) 

{  puts("***  Error  reading  data  file  ***"); 
exit (0)  ; 

} 

> 

void  mainQ 

{ 

read_input_f ile() ; 
read_dat_f ile() ; 
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mod()  ; 

> 


B.2  Shell  Programme 


'/,  This  section  converts  the  formatted  ASCI  files  into  a  binary  form 
for  use  by  SarTool 


cnvunf 

cnvunf 

cnvunf 

cnvunf 

cnvunf 

cnvunf 

cnvunf 

cnvunf 

cnvunf 

cnvunf 


test.fmt  0 

test.rfl . 0 

test.fmt . 1 

test.rfl . 1 

test.fmt . 2 

test.rfl . 2 

test.fmt .3 

test.rfl . 3 

test.fmt .4 

test.rfl .4 

test.fmt . 5 

test.rfl . 5 

test.fmt . 6 

test.rfl . 6 

test.fmt .7 

test.rfl .7 

test.fmt . 8 

test.rfl .8 

test.fmt . 9 

te'it.rfl .  9 

'/.Each  of  the  following  units  sequentially  copies  the  subject 
test_rfl  files  into  a  test.rfl  which  is  called  by  the  SarTool.dat 
file  and  used  as  the  subject  file  for  each  SarTool  run.  The 
resulting  output  is  placed  into  a  numbered  rpdb  file  and  then 
"stripped"  by  the  strip  programme.  The  output  of  the  strip 
programme  is  then  copied  to  a  test.out  file  for  later  use  by  the 
Mathematica  routine. 


cp  test.rfl. 0  test.rfl 
sartool  sartool.dat  rpdb.O 
dbstrip  rpdb.O  test.out. 0 

cp  test.rfl. 1  test.rfl 
sartool  sartool.dat  rpdb.l 
dbstrip  rpdb.l  test.out.l 

cp  test.rfl. 2  test.rfl 
sartool  sartool.dat  rpdb. 2 
dbstrip  rpdb. 2  test_out.2 
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cp  test_rfl.3  test.rfl 
sartool  sartool.dat  rpdb.3 
dbstrip  rpdb.3  test_out.3 

cp  test.rfl. 4  test.rfl 
sartool  sartool.dat  rpdb.4 
dbstrip  rpdb.4  test_out.4 

cp  test.rfl. 5  test.rfl 
sartool  sartool.dat  rpdb.5 
dbstrip  rpdb.5  test_out.5 

cp  test.rfl. 6  test.rfl 
sartool  sartool.dat  rpdb.C 
dbstrip  rpdb.6  test.cut.6 

cp  test.rfl. 7  test.rfl 
sartool  sartool.dat  rpdb.7 
dbstrip  rpdb.7  test_out.7 

cp  test.rfl. 8  test.rfl 
sartool  sartool.dat  rpdb.8 
dbstrip  rpdb.8  test_out.8 

cp  test.rfl. 9  test.rfl 
sartool  sartool.dat  rpdb.9 
dbstrip  rpdb.9  test_out.9 

/  The  final  line  of  the  shell  programme  removes  the  test.rfl  file 
in  preparation  for  the  next  SarTool  run. 

rm  test. .  rf  1 


B.3  Strip  Software 


The  strip  software  ’’strips”  the  output  of  SarTool  and  converts  it  into  in  phase 
and  quadrature  components  for  later  use  by  the  Mathematica  routine. 
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- - 

c 

c  program  dbstrip 

c 

c - 

c 

c  Routine  --  dbstrip  (standalone) 

c 

c 

c  Purpose  —  To  extract  the  magnitude  of  the  4  polarizations 

c  of  a  SarTool  range  profile  database  file 

c 
c 

c  Usage:  dbstrip  SarTool_database_f ile,  Output_file 

c 

c 

c  Command-line  arguments: 

c 

c 

c  (1)  Input  file  name  (example:  rpdb_mm_dd .xyz) 
c  implicit  undefined  (a-z) 

c  -  The  following  are  declarations  for  the  database  history 

variables 

double  precision  sldsvr,  rarexv,  toolv,  xsnsrO,  ysnsrO, 

zsnsrO,  fltang,  divang,  snsspd,  xtgt,ytgt, 
ztgt ,phit .thetat ,  psit,  vxtgt ,  vytgt , 
vztgt ,  axtgt ,  aytgt.aztgt,  freq,  dbwel, 
dbwaz,  dborel,  dboraz,  aspinc,  dr,  da, 
surfac,  bckatt,  thntrs,  dbter(4) ,  rtod 

integer  jdate(9),  irxdat(9),  nr,  na,  np,  i,  nintrv, 

imagtr,  iphstr  iargc,  nargs 

logical  isrflg,  sptflg,  shdflg,  terflg,  bshflg, 

cmpflg,  thnflg, 
wgtflg,  mltflg,  movflg 
character* 50  title 
character*30  rrxttl 
character*14  tgtnam 
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real*8  Ill,  Qll,  112,  Q12,  121,  Q21,  122,  Q22 
real*8  PAmp,  PPhase 

double  precision  Pll(500,2),  P12(500,2),  P21(500,2), 

P22 (500 , 2) 

integer  inputunit,  outputunit 
logical  create 
character*50  infil,  outfil 

c  -  Begin  executables 

data  rtod  /57 . 29577951/ 

c  -  Check  the  command  line  for  the  proper  number  of  arguments 

nargs  =  iargcQ 
if  (nargs  .It.  1)  then 
write(*  ,*) 

'  Usage:  dbstrip  SarTool_database_f ile , 
Output_text_f ile ' 
stop 
endif 

c  -  If  all  okay,  retrieve  the  command  line  argument 

call  getarg(l, infil) 
call  getarg(2, outfil) 

c  -  Define  the  SarTool  database  input  file  logical  unit  number 

inputunit  *  9 
outputunit  =11 

c  -  Open  the  file  containing  the  range  profile  database 

inquire  (file  =  infil,  exist  =  create) 
if  (.not.  create)  then 

write  (6,*) 'Input  file' .infil, 'does  not  exist.' 
stop 
endif 
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open  (unit  =  inputunit , 

file  =  infil, 

access  =  'direct  ’ , 
form  =  'unformatted', 
status  =  ' old' , 
reel  =  64) 

c  -  Read  parameters  from  the  database  header 


read ( inputunit , rec*  1 ) 

title. 

tgtnam 

read ( inputunit ,rec=  2) 

jdate 

read (inputunit ,rec=  3) 

rrxttl 

read ( inputunit ,rec=  4) 

sldsvr , 

rarexv .toolv 

read ( inputunit , rec=  5 ) 

irxdat 

read ( inputunit ,rec=  6) 

xsnsrO , 

ysnsrO , 

zsnsrO , 

f ltang , 

divang, 

snsspd, 

xtgt, 

Ytgt 

read ( inputunit ,rec=  7) 

ztgt. 

phit , 

thetat , 

psit , 

vxtgt , 

vytgt. 

vztgt , 

axtgt 

read ( inputunit ,rec=  8) 

aytgt , 

aztgt , 

freq. 

dbvel , 

dbwaz, 

dborel , 

dboraz. 

aspinc 

read ( inputunit ,rec=  9) 

dr. 

da. 

surf ac , 

bekatt , 

dbter(l),  dbter(2) , 

dbter(3) , 

dbter (4) 

read (inputunit ,rec=10) 

thntrs 

read (inputunit ,rec=ll) 

nr, 

na, 

np. 

nintrv, 

imagtr. 

iphstr, 

isrflg. 

sptflg. 

shdf lg, 

terflg. 

bshf lg. 

cmpflg. 

thnflg, 

wgtflg, 

mltf lg, 

movf lg 

DO  100,  i=12,  np+11 

read(inputunit ,  rec; 

=i)  Ill, 

Qll,  112,  Q12,  121,  Q21,  122, 

Q22 

Pll (i-11 , 1)=I11 
Pll(i-11,2)=Q11 
P12(i-ll , 1)=I12 
P12(i-11,2)=Q12 
P21(i-ll,l)=I21 
P21(i-11,2)=Q21 
P22(i-ii , 1)=122 
P22(i-11,2)=Q22 
100  continue 

close  (inputunit) 
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open  (outputunit,  file=outfil) 
do  200,  i=l,  np 


PAmp=(Pll(i,l)**2  +  Pll(i,2)**2)**.5 
IF  (Pll (i , 1) . EQ .0)  THEN 
PPhase=0 
ELSE 

PPhase=DATAN2(Pll(i , 2) , PI  1 (i , 1) ) 
END  IF 

write (outputunit ,900)  PAmp , PPhase 

PAmp=(P12(i,l)**2  +  P12(i ,2)**2)**.5 
IF  (P 12 (i , 1) . EQ . 0 )  THEN 
PPhase=0 
ELSE 

PPhase=DATAN2(P12(i,2) ,P12(i,l)) 
END  IF 

write (outputunit ,900)  PAmp .PPhase 

PAmp=(P21 (i , 1)**2  +  P21 (i ,2)**2)** . 5 
IF  (P21(i , 1) .EQ .0)  THEN 
PPhase=0 
ELSE 

PPhase=DATAN2 (P21 ( i , 2) , P2 1 ( i , 1) ) 
END  IF 

write (outputunit ,900)  PAmp, PPhase 

PAmp=(P22(i, 1)**2  +  P22(i,2)**2)**  .5 
IF  (P22(i , 1) .EQ. 0)  THEN 
PPhase=0 
ELSE 

PPhase=DATAN2(P22(i ,2) ,P22(i,l)) 
END  IF 

write(outputunit , 900)  PAmp , PPhase 

200  continue 

close  (outputunit) 

stop 

c  -  Format  statements 
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900  f ormat (el2 . 5) 


end 


B.4  Mathtmaiitu  Rcutine 


The  following  programme  is  an  abridged  version  of  the  Mathematica  routine. 
Each  intermediate  unit  of  the  routine  is  preceeded  by  functional  documentation. 


The  following  two  files  contain  the  "source  code"  for  the 
parametric  plots  and  the  Poincare  sphere 

The  next  line  reads  in  the  In  Phase  and  Quadrature  data  from 
SarTool  using  the  "  Strip  Programme". 

fulllist=ReadList  [Inputstring  ["input  file  name"] , 

{{Real,  Real},  {Real,  Real},  {Real,  Real},  {Real,  Real}}]; 

interval=7 ; 
np=200 ; 

pulses=Floor [np/interval] ; 
m={}; 

Block[  {i}, 

m=Table[fulllist [[i] ] ,{i,l,np+l-interval, interval}] ; 

]; 

The  following  four  tables  tabulate  the  Scattering  Matrix 

components 

for  each  pulse 

saa=Table[(m[[i, 1 ,l]]*Exp[I  m[[i, 1 ,2]]] ) ,{i, pulses}] ; 
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sba=Table[(m[[i,3,l]]*Exp[I  m[[i,3,2]]]) ,{i .pulses}] ; 
sbb=Table[(m[[i,4,l]]*Exp[I  m[[i,4,2]]]) ,{i .pulses}] ; 
sab=Table [(m[[i ,2 , 1] ] *Exp[I  m[[i,2,2]]]) ,{i .pulses}] ; 


The  next  two  lines  calculate  the  intermediate  variables  rho  1(22 
(pic  <8  p2c)  which  are  used  to  calculate  qlc  4  q2c  which  are  in  turn 
used  for  the  eventual  calculation  of  the  latitude  &  longitude 
of  the  copols  on  the  Poincare  Sphere. 

plc=  Tablet  If [sbb [ [i] ] ==0 , 

If [sab [ [i] ] ==0,1,1] , 

(-(sab[[i]]) 

+  Sqrt [(sab [ [i]] ) '2 

-  sbb[[i]]*saa[[i]]])/ (sbb [ [i] ] ) ] 

,{i, pulses}] ; 

p2c=  Tablet  If [sbb[[i]]==0, 

If [sab [ Ci] ] ==0,-1 , -Sign [sab [ [i] ]  ] 
♦Infinity] , 

(-(sab[[i]]) 

-  Sqrt[(sab[[i]])~2 

-  sbb[[i]]*saa[[i]]])/(sbb[[i]])] 

, {i  .pulses}] ; 


The  following  section  determines  the  values  of  rl,r2  ®  r3  which 
are  intermediate  results  used  to  determine  the  variables  rho  102 
(plx  0  p2x)  which  in  turn  are  used  to  calculate  the  latitude  and 
longitude  of  the  crosspols  on  the  Poincare  Sphere. 

rl=  Tablet  (  (Abs [sbb[[i]]] ) ‘2 

-(Abs[saa[[i]]])“2) 

-sab [ [i] ] *Conj ugate [sba [ [i] ] ] 

+sba [ [i] ] *Conjugate [sab  [  [i] ] ] 

, {i, pulses}] ; 


r2=  Tablet  (  sbb[[i]] ♦Conjugate[sba[[i]]] 

+sba [ [i] ] *Conj ugate [saa [ [i] ] ] 

, {i, pulses}] ; 
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r3=  Tablet  (  saa[ [i] ] *Conjugate [sab [ [i]]] 

+sab[[i]] *Conj  ugate [sbb  [  [i] ] ] 


,{i, pulses}] ; 


) 


plx=  Table[  If [r2[[i]]==0, 

if [rl [ti]]>01Sign[rl[[i]]]*Infinity,0] , 

(  rl[[i]] 

+  Sqrt [r 1  [  [i J ] *2 

+  4*r2[[i]]*r3[[i]]])/(2*r2[[i]])] 

, {i, pulses}] ; 

p2x=  Tablet  If [r2 t [i] ] ==0 , 

If [rl tti]] <0,Sign[rl[[i]]]*Inf inity ,0]  , 

(  rlCCi]] 

-  Sqrt trl [ti]] “2 

+  4*r2 [ [i] ] *r3  [[i] ] ] )/(2*r2[[i]])] 

, {i .pulses}] ; 


The  following  sections  determine  the  values  of  q  for  each  of  the 
p  values.  These  q  values  are  subsequently  used  to  determine  the 
theta  and  phi  values  of  the  cross  and  copol  nulls  on  the 
Poincare  Sphere. 

qlx=  Table [If  [SameQ [plx tti]] .Infinity]  II  SameQ tplx t ti] ] , 

-Infinity] 

.-1 

.If [SameQ[plx[[i]] ,1.  I] 

,1 

, (1-I*plxt[i]l )/ (1+I*plx[[i]] )]] 

,{i, pulses)]  ; 

q2x=TabletIf  [SameQ [p2x tti]] .Infinity]  If  SameQ [p2x[[i]]  , 

-Infinity] 

,-l 

,If tSameq[p2xtti]] ,1.  I] 

.1 

, (1-I*p2x[[i]] )/(l+I*p2xtti]])]] 

,{i, pulses}] ; 

qlc=Table [If [SameQ tplct ti]] , 1 .  I] 

,1 
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,((1-I*plc[[i]])/(1+I*plc[[i]]))] 
,  {i .pulses)] ; 


q2c=Table [If  [SameQ[p2c [  [i] ] , 1 .  I] 

,1 

, ( (1-I*p2c[  [i] ] )/ (1+I*p2c  [[i]] ) )] 

,{i .pulses}]  ; 

The  following  section  determines  the  co-latitude  for  each 
of  the  above  q  values. 

crosspolelLat=Table [N[ArcCos [- ( ( 1 - ( Abs [qlx  [[i]]])“2) 

/(l+(Abs[qlx[[i]]])~2))] ,5] 

,{i .pulses}] ; 

crosspole2Lat=Table [N[ArcCos [-( (l-(Abs [q2x [ [ i] ] ] ) ‘2) 

/ (1+ (Abs [q2x[ [i] ]])*2))] ,5] 

,{i, pulses}]  ; 

copolelLat*Table[If  [SameQ  [plc[[i]]  ,1.1] 

,0 

,N[ArcCos[-((l-  bs [qlc[[i] ]] ) ~2) 
/(l+(Abs[qlc[[i]]])-2))]],5] 

,{i .pulses}]  ; 

copole2Lat=Table [If [SameQ [p2c [ [i] 3 , i . I] 

,0 

,N [ArcCos [- ( ( 1- (Abs [q2c [ [i] ] ] ) *2) 
/ (l+(Abs[q2c[[i]]] )*2))]] ,5] 

,{i, pulses}] ; 


The  following  section  determines  the  longitude  of  each  of  the  cross 
and  copol  nulls. 

crosspolelLong=Table[lf [Re [qlx [[i]]] ==0 , 

If  [lm[qlx[[i]]]«0, 

N[ArcTan[-l  .0]  ,5]  , 
-Sign[Im[qlx[[i]]]]*N[Pi/2,5]] , 
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,{i .pulses}] ; 


N [ArcTan [Re [qlx [ [1] ] ] ,-Im[qlx[[i]]]] 
]  .5  ] 


crosspole2Long=Table [If [Re[q2x [[!]]] ==0, 

If [Im[q2x[[i]]]==0, 

N[ArcTan[-l  .0]  ,5]  , 

-Sign[Im [q2x [ [i]] ] ] *N  [Pi/2 , 5] ] , 

N [ArcTan [Re [q2x [ [i] ] ] , -Im [q2x  [  [l] ] ] ] 
]  ,  5  ] 

,{i .pulses}] ; 

copole lLong=Table [If [Re [qlc [ [i] ]] ==0 , 

If [lm[qlc[[i]]]==0, 

N [ArcTan [-1 . 0] , 5] , 

-Sign[Im [qlc [ [i]] ] ] *N [Pi/2 , 5] ] , 

N [ArcTan [Re [qlc [ [i] ] ] . -Im [qlc [ [i] ] ] ] 
],5] 

,{i .pulses}] ; 

copole2Long=Table [If [Re [q2c [ [i] ] ] ==0 , 

If [Im[q2c[[i]]]=*0, 

N [ArcTan [- 1 . 0] , 5]  , 

-Sign[Im[q2c [[i]] ] ] *N[Pi/2 ,5] ] , 

N [ArcTan [Re [q2c [ [i] ] ] , -Im [q2c [ [l] ] ] ] 
],5  ] 

,{i, pulses}] ; 


The  following  section  plots  the  co  and  cross  polarizations  on  the 
Poincare  Sphere. 

copoleList=Table[{copolelLat [[i]] ,  copole 1 Long [ [i] ] } ,{i .pulses}] ; 

Do [AppendTo [copoleList , {copole2Lat [ [i] ] , 
copole2Long[[i]] }] 

, {i .pulses}] ; 

crosspoleList=Table [{crosspolelLat [ [i] ] , 
crosspole lLong[[i]]} ,{i .pulses}] ; 
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Do  [AppendTo  [crosspoleList  , {crosspoie2Lat  [  [i]  ] 

,  crosspole2Long[[i]]}] 

,{i ,pulses>]  ; 

The  following  is  appended  to  the  Mathematica  sphere  graphic 

Thickness [.009] ,  Line[{{0. ,0. , l.},{0. ,0. ,-!.}>], 

Thickness [. 005] ,  Line [{{1. ,0, ,0. >,{-1 . ,0. ,0 . }}] , 

Thickness [ . 001] ,  GrayLevel  [.5] .Line [{{0. ,1. ,0.}, 

{0  .  ,-l.  ,0.»]  ,  GrayLevel[0]  >; 
copolelen  =  Length[copoleList] ; 
copoleplot=Table [{PointSize [ . 02] , 

Point [{Sin[copoleList  [ [i , 1] ] ] 

Cos [copoleList  [[i ,2]]]  , 

crosspolelen  =  Length [crosspoleList] ; 
crosspoleplot=Table[  {PointSize[.02]  , 

GrayLevel [. 5] .Point [{Sin [crosspoleList  [ [i  ,  l] ] ] 
Cos[crosspoleList[[i,2]]] , 

AppendTo [sphere .copoleplot] ; 

AppendTo [sphere, crosspoleplot] ; 

The  following  graphics  shows  the  co  and  cross  polarizations  plotted 
on  the  Poincare  sphere.  The  copolarizations  are  the  darker  of  the 
points.  Refer  to  main  body  of  thesis  for  examples  of  this  output. 

test=Graphics3D[sphere] ; 

Show[test,  Boxed  ->  False]; 

The  following  graph  shows  the  span  of  the  returns. 

span  =  Table[(Abs[saa[[i]]] “2  +  Abs[sba[[i]]]*2  +Abs [sab  [[i] ] ] '2 

+  Abs[sbb[[i]]]~2),{i,l, pulses)] 

ListPlot [crosspolelLat .Plot Joined  ->  True,  PlotRange  ->{-3.2, 3. 2}, 

PlotLabel  ->  "crosspolel  Lat"] 

testl=ListPlot [crosspolelLong, Plot Jo ined->True .PlotRange- > {-3 .2,3 
•  2}, 


B-38 


PlotLabel  -V'crosspolel  Long"] 


List Plot [copole lLat , Plot Joined->True ,PlotRange->{-3 . 2 ,3 . 2}  , 

PlotLabel  ->  "copolel  Lat"] 

ListPlot [copol el Long .Plot Jo ined->True, Plot Range -> {-3 .2,3.2}, 

PlotLabel  ->"copolel  Long"] 

ListPlot [crosspole2Lat .Plot Joined  ->  True,  PlotRange  ->{-3 . 2 ,3 . 2} , 

PlotLabel  ->  "crosspole2  Lat"] 

test2=ListPlot [crosspole2Long , Plot Joined- >True .PlotRange- > {-3 .2,3 

•  2} , 

PlotLabel  ->"crosspole2  Long"] 

ListPlot [copole2Lat .Plot Joined- >True ,PlotRange->{-3. 2,3 . 2} , 

PlotLabel  ->  "copole2  Lat"] 

ListPlot [copole2Long,Plot Joined->True ,PlotRange->{-3 .2 ,3 . 2} , 

PlotLabel  ->  "copole2  Long"] 

z=Table[{copolelLat[[i]] , copolelLong[[i]] } ,{i , 1 .pulses}] 

ListPlot [z,  PlotLabel->"copolelLat  versus  copolel . Long"] 
y=Table [{copole2Lat [ [i] ] ,copole2Long[[i]] } ,{i , 1 .pulses}] 
Show[testl,test2] 


B.5  Convert  Programmme 


The  following  programme  converts  the  mathematica  output  into  a  form  suitable 
for  input  to  G.  Tarr’s  multi  level  perceptron  network. 


DIM  elem$(l00) 

INPUT  "What  file  contains  the  filenames  to  be  processed"  ,  n 


OPEN  "n"  FOR  INPUT  AS  #1 
DO  UNTIL  EOF ( 1) 

INPUT  #1,  file$ 

OPEN  file$  FOR  BINARY  AS  #2 
OPEN  "dummy.dat"  FOR  BINARY  AS  #3 
DO  UNTIL  EOF (2) 

GET  #2,  elem 

IF  elemS  =  OR  elem$  =  THEN 
elemS  =  STRING$(1,  0) 
ELSEIF  elemS  =  THEN 

elemS  =  STRINGS (1,  13) 

PUT  #3,  elem 
elemS  =  STRINGS (1,  10) 
PUT  #3,  elem 

ELSE 

PUT  #3,  elem 

END  IF 
LOOP 

CLOSE  #2,  #3 
LOOP 

KILL  file! 

NAME  "dummy.dat"  AS  file$ 

END 


B.  6  Single  Perceptron  Network 

The  following  programme  creates  a  single  level  perceptron.  The  basis  of  this 
programme  was  written  for  an  AT.  unit  by  the  author  and  Capt  J.  De-Barry.  The 
author  then  modified  the  programme  such  that  it  could  train  and  test  on  opera¬ 
tor  selected  inputs  and  return  decision  files  and  outputs.  Extensive  documentation 
appears  at  the  end  of  the  programme. 
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’DECLARE  SUB  instruct  () 

'DECLARE  SUB  init  (eta!,  weight# (),  theta!,  count!,  right!,  wrong1, 
maxcount ! ,  0$,  s!,  train! ,  n!,  r,  outfile$,  trainfile$) 

'DECLARE  SUB  distl  (vector#(),  class!) 

'DECLARE  SUB  dist2  (vector#(),  class!) 

'DECLARE  SUB  dist3  (vector#(),  class!) 

'DECLARE  SUB  dist4  (vector#(),  class!) 

'DECLARE  SUB  readin  (DATAFILE!,  vect!(),  class !(),  n! ,  trainfileS) 
'DECLARE  SUB  dist5  (vector#(),  class!,  class !(),  vect!(),  count!, 
n !  ) 

'DECLARE  SUB  feedforward  (guess!,  vector#(),  weight#(),  theta!, 
count!,  train!,  maxcount!,  x!,  y2) 

'DECLARE  SUB  backprop  (guess!,  class!,  eta!,  weight#(),  vec tor#(), 
theta ! ,  y2) 

'DECLARE  SUB  backprop2  (guess!,  class!,  right!,  wrong!,  ratio!, 
weight#() ) 

'DECLARE  SUB  display  (vector#(),  ratio!,  weight#()  ,  theta!,  d!, 
guess!,  class!,  count!,  maxcount!,  s!,  train!,  r,  outfileS,  0$) 

DIM  vector#(l),  weight#(l) 

DIM  vect (450 ,  2) 

DIM  class(450) 

RANDOMIZE  TIMER 

DEF  FNyl#  (m,  x,  B)  =  (-m)  *  x  +  B 

SCREEN  9 
CALL  instruct 
DO 

CALL  init(o ,q,eta,  veight#(),  theta,  count,  right,  wrong, 
maxcount,  0$,  s,  train, n,r,p  ) 

WINDOW  (-6.5,  -6.5)-(6 .5,  6.5) 

VIEW  SCREEN  (PMAP(-3.25  +  r,  0),  PMAP(3.25  +  s , 1) ) - (PMAP(3.25 
+  r,  0)  ,  PMAP(-3.25  +  s,  D) 

blockstep  =  maxcount 
DO 

IF  0$  =  "1"  THEN 

CALL  distl(vector#() ,  class) 

ELSEIF  0$  =  "2"  THEN 

CALL  dist2(vector#() ,  class) 

ELSEIF  0$  -  "3"  THEN 


CALL  dist3(vector#() ,  class) 

ELSEIF  C$  =  "4"  THEN 

CALL  dist4(vector#() ,  class) 

ELSEIF  0$  =  "5"  THEN 

IF  count  =  1  THEN 

CALL  readin(q .DATAFILE,  vectQ,  classO,  n, 
trainf ile$ ,tsstfile$) 

ELSEIF  count  >  1  THEN 

CALL  dist5(k,vector#()  ,  class,  class(),  vectQ, 
count,  n,o) 

END  IF 

END  IF 

CALL  feedforward (guess ,  vector#()  ,  weight#(),  theta, 
count,  train,  maxcount,  x,  y2) 

IF  train  =  maxcount  OR  count  <  train  THEN 

CALL  backprop (guess ,  class,  eta,  veight#(), 
vector#(),  theta,  y2) 

ELSEIF  train  <>  maxcount  AND  count  >  train  THEN 
IF  count  =  train  +1  AND  qOl  AND  0$="5"  THEN 
CALL  readin2, DATAFILE,  vect (), class (), o , 
trainf ile$,testfile$) 

END  IF 

CALL  backprop2 (guess,  class,  right,  wrong,  ratio, 
weight# () ) 

END  IF 

CALL  display(k,p,vector#() ,  ratio,  weight#(),  theta,  d, 
guess,  class,  count,  maxcount,  s,  train,  r, 
outfile$,  0$) 
count  =  count  +  1 
IF  count  =  maxcount  +  1  THEN 
DO 

LOCATE  23,  8 

IF  train  *  maxcount  OR  train  >  maxcount  THEN 

INPUT  "Do  you  want  to  continue  training  (Y/N)" ;  A$ 

ELSEIF  train  <  maxcount  THEN 

INPUT  "Do  you  want  to  continue  testing  (Y/N)";  A$ 

END  IF 

LOOP  UNTIL  A$  =  "y"  OR  A$  =  "Y"  OR  A$  ■  "n"  OR  A$  =  "N" 
LOCATE  23,  6 
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PRINT 


IF  A$  =  "y"  OR  A$  =  "Y"  THEN 

maxcount  =  maxcount  +  blockstep 
LOCATE  1,  25 
PRINT  "DYNAMIC" ; 

END  IF 
END  IF 

LOOP  UNTIL  count  =  maxcount  +  1 
LOCATE  23,  8 

INPUT  "Do  you  want  run  the  simulation  again  (Y/N)";  A$ 
CLS 

WINDOW 

VIEW 

CLS 

LOOP  UNTIL  A$  =  "n"  OR  A$  =  "N" 


END 

SUB  backprop  (guess,  class,  eta,  weight#(l),  vector#(l),  theta,  y2) 

'  IF  guess  <>  class  THEN 

FOR  I  =  0  TO  1 

veight#(I)  =  weight#(I)  +(eta/2)  *  (class  -  y2)  * 

vector#(I) 

NEXT  I 

theta  =  theta  -  eta  /  2  *  (class  -  y2) 

'  END  IF 

'  LOCATE  13,1 

'  PRINT  "b/l  weight#(0)="  veight#(0) 

’  PRINT  "b/l  weight#(l)="  veight#(l) 

END  SUB 

SUB  backprop2  (guess,  class,  right,  wrong,  ratio,  weight#(l)) 

'  LOCATE  15,1 

’  PRINT  "b/2  weight#(0)=  "  weight#(0) 

’  PRINT  "b/2  weight#(l)=  "  weight#(l) 

IF  guess  <>  class  THEN 
wrong  *  wrong  +  1 
ELSEIF  guess  *  class  THEN 
right  *  right  +  1 
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END  IF 

ratio  =  right  /  (right  ♦  wrong) 

LOCATE  10,  1 
PRINT  "Testing  Data" 

PRINT  "wrong=" ;  wrong 
PRINT  "right=";  right 
PRINT  "'/.corrects";  ratio  *  100 
END  SUB 

SUB  display  (k,p, vector#(l) ,  ratio,  weight#(l),  theta,  d,  guess, 
class,  count,  maxcount ,  s,  train,  r,  outfile$,  0$) 

SHARED  trainfile$,  testfile$ 

STATIC  oldm,  oldb 
IF  count  =  1  THEN 
LOCATE  1,  25 

PRINT  "INITIAL  PERCEPTRON  DECISION  SPACE"; 

LOCATE  2,  23 
IF  0$  =  "5"  THEN 

PRINT  "INPUT  TRNG/TST  FILES  "; 
trainf ile$; "/" ;testf ile$ 

END  IF 
oldm  =  1 
oldb  =  .2 

ELSEIF  count  =  2  THEN 
LOCATE  1,  25 
PRINT  "DYNAMIC"; 

ELSEIF  count  +  1  >=  maxcount  THEN 
LOCATE  1,  25 
PRINT  "  FINAL"; 

END  IF 

IF  count  <  train  THEN 
IF  0$*"5"  THEN 
LOCATE  3,  34 
PRINT  "TRAINING  PHASE" 

ELSEIF  0$<>"5"  THEN 
LOCATE  2,34 
PRINT  "TRAINING  PHASE" 

END  IF 

ELSEIF  count  >  train  THEN 
IF  0$*"5"  THEN 
LOCATE  3,27 

PRINT  "TESTING  PHASE  -  WEIGHTS  FROZEN" 
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ELSEIF  0$<>"5"  THEN 
LOCATE  2,27 

PRINT  "TESTING  PHASE-  WEIGHTS  FROZEN" 

END  IF 
END  IF 
LOCATE  4,  1 

PRINT  "vector#(0)=" ;  CSNG(vector#(0)) 

PRINT  "vector#(l)=" ;  CSNG(vector#(l) ) 

PRINT  "class=";  class 
PRINT  "guess=";  guess 

LINE  (-3.2  +  r,  FNyl#(oldm,  -3.2,  oldb)  +  s)-(3.2  + 
r ,FNyl#(oldm,  3.2,  oldb)  +  s) ,  0 

LINE  (0  +  r,  3.2  +  s)-(0  +  r,  -3.2  +  s) ,  15 
LINE  (-3.2  +  r,  0  +  s)-(3.2  +  r,  0  +  s),  15 
newm#  =  weight#(0)  /weight#(l) 
newb#  =  theta  /veight#(l) 

LINE  (-3.2  +  r,  FNyl#(newm#,  -3.2,  newb#)  +  s)-(3.2  + 
r ,FNyl#(newm#,  3.2,  newb#)  +  s) ,  2 
IF  count  =  1  THEN 
LOCATE  5,  46 
PRINT  "3.2"; 

LOCATE  17,  45 
PRINT  "-3.2"; 

LOCATE  11,  25 
PRINT  "-3.2"; 

LOCATE  11,  68 
PRINT  "3.2"; 

LOCATE  23,  8 

INPUT  "Press  enter  to  begin  training  A$ 

LOCATE  23,  8 
PRINT  " 

END  IF 

IF  class  *  1  THEN 

PSET  (vector#(0)  +  r,  vector#(l)  +  s) ,  14 
ELSEIF  class  »  -1  THEN 

PSET  (vector#(0)  +  r,  vector#(l)  +  s) ,  9 
END  IF 

oldm  *  weight#(0)  /  weight#(l) 
oldb  *  theta  /  weight#(l) 

LOCATE  19,  1 

PRINT  ”W(0)  =  ";  CSNG (weight# (0) ) ; 

LOCATE  19,  40 
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PRINT  "W(l)  =  CSNG (weight# (1)) 

LOCATE  20,  20 

PRINT  "Theta  =  " ;  CSNG(theta) 

LOCATE  21,  1 

PRINT  "ITERATION  -  ";  count; 

LOCATE  21,  40 

PRINT  "INTERCEPT  =  ";  CSNG(oldb) 

LOCATE  22,  20 

PRINT  "SLOPE  =  ";  CSNG(oldm) 

LOCATE  23,  8 
PRINT  " 

corr#  =  ratio*100 

IF  0$  =  "5"  AND  p»l  AND  count  >  train  THEN 
OPEN  outfile$  FOR  APPEND  AS  #2 
IF  count  =  train  +  1  THEN 

PRINT  #2,  "  " ; outf ile$ ; " 


PRINT  #2, 

PRINT  #2,"  Selection  Class  Guess 

'/.Correct" 

PRINT  #2,  SPC(9)  count  SPC(9)  class  SPC(9)  guess  SPC(9) 
corr# 

ELSEIF  count  >  train  +  1  THEN 

PRINT  #2,  SPC(9)  count  SPC(9)  class  SPC(9)  guess  SPC(9) 
corr# 

END  IF 
CLOSE  #2 
END  IF 


END  SUB 

SUB  distl  (vector#(l),  class) 

vector#(0)  =3*(RND  *2-1) 
IF  vector# (0)  >  0  THEN 
vector#(l)  *  RND  *3 
class  *  1 

ELSEIF  vector#(0)  <  0  THEN 
vector# (1)  =  RND  *  (-3) 
class  =  -1 
END  IF 
END  SUB 
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SUB  dist2  (vector#(l),  class) 

IF  RND  >  .5  THEN 
class  =  1 

vector#(0)  *  3*(  RND  *  (-.9)  -  -1) 
vector#(l)  =  3*(RND  *  (.6)  +  .1  ) 

ELSE 

class  =  -1 

vector# (0)  =3*(RND  *  1.2  -  .6) 

IF  vector# (0)  <  .3  THEN 

vector#(l)  =  3* (RND  *  (-.4)  -  .1) 
ELSE 

vector#(l)  =3* (  RND  *  1.2  -  .5) 

END  IF 

END  IF 
END  SUB 

SUB  dist3  (vector#(l),  class) 
vector#(0)  =  RND  *  6  -  3 
IF  vector# (0)  >  .6  THEN 

vector#(l)  *  3* (RND  -  .2) 
class  *  1 

ELSEIF  vector#(0)  <  -.6  THEN 

vector#(l)  =  3* (RND  *  (-1)  +  .2) 
class  =  -1 
ELSE 

IF  RND  >  .5  THEN 

vector# (1)  *  3* (RND  -  .2) 
class  *  1 
ELSE 

vector# (1)  =  3* (RND  *  (-1)  +  .2) 
class  *  -1 
END  IF 
END  IF 
END  SUB 

SUB  dist4  (vector#(l),  class) 

IF  RND  >  .5  THEN 
class  a  1 

vector#(0)  »3*(RND  *  (-.9)  -  .1) 
vector#(l)  *3*(RND  *  (.6)  +  .1) 
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ELSE 

class  =  -1 

vector#(0)  =3*(RND  -  .4) 

IF  vector# (0)  <  .3  THEN 

vector#(l)  *  3*(RND  *  (--4)  -  .1) 

ELSE 

vector#(l)  =  3* (RND  *  .8  ~  .7) 

END  IF 

END  IF 
END  SUB 

SUB  dist5  (k,vector#(l) ,  class,  class(l),  vect(2),  count,  n,o) 
IF  trainf ile$  =  testfilel  THEN 
DO  UNTIL  k  <>  0 
k  =  CINT(RND  *  n) 

LOOP 

ELSEIF  trainf  ile$Otestfile$  THEN 
DO  UNTIL  k<>0 
k  =  CINT(RND  *  o) 

LOOP 
END  IF 
LOCATE  8,1 

PRINT  "vector  selected="  k 
vector#(0)  =  vect(k,  1) 
vector#(l)  =  vect (k,  2) 
class  =  class (k) 
k  =  0 

END  SUB 

SUB  feedforward  (guess,  vector#(l),  veight#(l),  theta,  count, 
train,  maxcount,  x,  y2) 

Z  *  vectort(O)  *  weight#(0)  +  vector#(l)  *  weight#(l) 

IF  Z  <  -1  OR  Z  >  1  THEN 
IF  Z  >  0  THEN 
guess  *  1 
y2  »  1 

ELSEIF  Z  <  0  THEN 
guess  *  -1 
y2  »  -1 


theta 
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END  IF 

ELSEIF  Z>0  THEN 
guess  =1 
y2  =Z 

ELSEIF  Z<0  THEN 
guess  =  -1 
y2  =  Z 

END  IF 

LOCATE  15,  1 

PRINT  "Initialization  params" 

PRINT  "train3" ;  train 
PRINT  "total  iterations3" ;  maxcount 
END  SUB 

SUB  init  (o,q,eta,  weight#(l),  theta,  count,  right,  wrong, 
maxcount,  0$,  s,  train, n,r,p  ) 

SHARED  outfile$,  trainfile$,  testfileS 

r  =  1! 

s  =  1 ! 

count  =  1 

wrong  *  0 

right  *  0 

weight#(0)  =  RND  *  .2 
weight#(l)  *  (RND  -  1)  *  .2 
theta  =  (RND  -  1)  *  .1 
DO 

LOCATE  10,  18 

PRINT  "Choose  class  separation  type  1,  2,  3,  4  or 

LOCATE  12,  21 

PRINT  "(1)  linearly  separable  (quadrant  I  and  III) 

LOCATE  13,  21 

PRINT  "(2)  linearly  unseparable" ; 

LOCATE  14,  21 

PRINT  "(3)  overlapping"; 

LOCATE  15,  21 

PRINT  "(4)  linearly  separable  (arbitrary)"; 

LOCATE  16,  21 

PRINT  "(5)  read  in  data"; 

LOCATE  18,  18 

INPUT  "Your  choice";  0$ 


LOOP  UNTIL  0$  =  "1"  OR  0$  =  "2"  OR  0$  =  "3"  OR  0$  =  "4"  OR 
0$  =  "5" 

CLS 

LOCATE  12,  24 

INPUT  "How  many  total  iterations";  maxcount 
CLS 

LOCATE  12,  14 

INPUT  "Do  you  wish  to  test  as  well  as  train  the  network 
(Y/N)";  A$ 

CLS 

IF  A$  =  "y"  OR  A$  =  "Y"  THEN 
LOCATE  12,24 

INPUT  "How  many  training  iterations";  train 
CLS 

ELSEIF  A$  =  "n"  OR  A$  =  "N"  THEN 
train  =  maxcount 
END  IF 

IF  0$  =  "5"  THEN 
LOCATE  12,  24 

INPUT  "What  training  input  file  name";  trainfileS 
CLS 

LOCATE  12,24 

INPUT  "What  testing  input  file  name";  testfile$ 

IF  trainfile$  =  testfileS  THEN 
q=l 

END  IF 
CLS 

LOCATE  12,  24 

INPUT  "How  many  entries  on  training  file";  n 
CLS 

IF  q<>l  THEN 
LOCATE  12,24 

INPUT  "  How  many  entries  on  testing  file";o 
CLS 

END  IF 

LOCATE  12,  18 

INPUT  "Do  you  wish  to  write  an  output  file  (Y/N)";P$ 
CLS 

IF  P$  =  "n"  or  P$  =  "N"  THEN 

p=0 

ELSEIF  P$  -  "y"  or  P$  »  "Y"  THEN 


P=1 


LOCATE  12,24 

INPUT  "What  output  file" 

CLS 

END  IF 

ELSEIF  0$  <>  "5"  THEN 
LOCATE  12,  24 
PRINT  SPC(65) ; 

END  IF 

LOCATE  12,  17 


outf ile$ 


INPUT  "What  training  gain  value  (0.3  suggested)"; 
CLS 

END  SUB 


eta 


SUB  instruct 
CLS 

LOCATE  8,  24 

PRINT  "  PERCEPT  " ; 

LOCATE  10,  23 

PRINT  "THE  PERCEPTRON  TRAINING  SIMULATOR" ; 

LOCATE  12,  35 
FOR  I  =  1  TO  6000 
NEXT  I 
CLS 
DO 

LOCATE  12,  25 

INPUT  "Do  you  want  instructions";  r$ 

LOOP  UNTIL  r$  =  "y"  OR  r$  *  "Y"  OR  r$  *  "n"  OR  r$  =  "N" 

CLS 

IF  r$  =  "y"  OR  r$  =  "Y"  THEN 
LOCATE  2,  4 

PRINT  "TRAINING  PARAMETERS" 

LOCATE  4,  6 

PRINT  "You  will  be  prompted  to  enter  several  training 
parameters . " 

LOCATE  5,  6 

PRINT  "First  you  will  be  asked  to  select  the  type  of  class 
separation" 

LOCATE  6,  6 

PRINT  "you  want.  Next  you  will  be  prompted  for  the  number 
of  total" 

LOCATE  7,  6 
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PRINT  "iterations  to  run  before  pausing.  You  may  input  any 
positive" 

LOCATE  8,  6 

PRINT  "integer.  You  will  then  be  asked  whether  or  not  you 
wish  to" 

LOCATE  9,  6 

PRINT  "perform  testing  with  the  net.  A  positive  response 
to  this" 

LOCATE  10,  6 

PRINT  "input  will  cause  the  weights  to  be  frozen  at  the 
specified" 

LOCATE  11,  6 

PRINT  "training  number  after  which  the  net  will  use  those 
weights  to" 

LOCATE  12,  6 

PRINT  "make  decisions  on  subsequent  inputs.  If  you 
initially  " 

LOCATE  13,  6 

PRINT  "selected  option  5  (read  in  data)  you  will  be  prompted 
to  input" 

LOCATE  14,  6 

PRINT  "the  number  of  entries  to  be  read  from  the  input  file. 
You  will" 

LOCATE  15,6 

PRINT  "also  be  asked  to  enter  the  input  and  output  file 
names.  With" 

LOCATE  16,6 

PRINT  "selection  5  there  is  also  the  capability  to  test  from 
a  file 
LOCATE  17,6 

PRINT  "other  than  the  one  trained  on  this  is  achieved  simply 
by 

LOCATE  18,6 

PRINT  "answering  the  questions  given  during  program 
initialization  Finally," 

LOCATE  23,  6 

INPUT  "Press  <RETURN>  to  continue  A$ 

CLS 

LOCATE  2,  6 

PRINT  "You  will  be  prompted  for  a  training  gain.  The 
suggested  value  " 
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LOCATE  3,  6 

PRINT  "is  0.3,  any  number  between  0.5  and  1  appears  to  work. 
Larger" 

LOCATE  4,  6 

PRINT  "values  train  faster,  but  the  decision  boundary  may 
oscillate" 

LOCATE  5,  6 

PRINT  "around  an  optimal  solution  when  input  classes  overlap. 

Small  " 

LOCATE  5,  6 

PRINT  "values  train  slower,  but  the  decision  boundary  will 
transition" 

LOCATE  7,  6 

PRINT  "smoothly  to  a  fairly  stable  solution.  The  program  and 
output" 

LOCATE  8,  6 

PRINT  "have  been  structured  such  that  the  operation  and 
results  are" 

LOCATE  9,  6 

PRINT  "readily  apparent.  For  example,  at  the  end  of  the 
training” 

LOCATE  10,  6 

PRINT  "phase  the  operator  is  alerted  to  the  fact  that  the 
weights" 

LOCATE  11,  6 

PRINT  "have  been  frozen  and  that  the  net  is  operating  in  the 
test" 

LOCATE  12,  6 

PRINT  "mode.  All  relevant  data  is  then  displayed  and 
updated  " 

LOCATE  13,  6 

PRINT  "including  the  ratio  of  correct  decisions  made  by  the 
net" 

LOCATE  15,  6 

PRINT  "The  following  page  describes  some  of  the  options 
available" 

LOCATE  16,  6 

PRINT  "to  the  operator" 

LOCATE  23,  6 

INPUT  "Press  <RETURN>  to  continue  ...  ",  A$ 

CLS 
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LOCATE  2,  4 
PRINT  "OPTIONS" 

LOCATE  4,  6 

PRINT  "There  are  several  options  available.  As  previously 
mentioned, " 

LOCATE  5,  6 

PRINT  "you  may  opt  to  continue  training  or  testing  after  the 
specified" 

LOCATE  6,  6 

PRINT  "number  of  iterations  has  been  reached.  If  you  opt  to 
continue , " 

LOCATE  7,  6 

PRINT  "the  net  will  start  at  the  current  iteration  count  and 
continue  " 

LOCATE  8,  6 

PRINT  "until  another  set  of  the  specified  iterations  is 
complete.  You  may" 

LOCATE  9,  6 

PRINT  "again  opt  to  continue  operation.  Once  you  quit 
operating,  you  may" 

LOCATE  10,  6 

PRINT  "opt  to  re  initializw  and  train  or  test  another  net 
under  different" 

LOCATE  11,  6 

PRINT  "parameters , if  you  choose  not  to  train  another  net,  the 
program  terminates." 

LOCATE  12,  6 

PRINT  "One  other  useful  option  is  provided.  This  program  is 
meant  to" 

LOCATE  13,  6 

PRINT  "stand  alone  and  does  not  contain  any  printer  routines. 
However,  if" 

LOCATE  14,  6 

PRINT  "you  have  the  EGA  screen  dump  utility  ’EGADMP',  you  can 
load  it  into" 

LOCATE  15,  6 

PRINT  "memory  first  (it  is  a  TSR)  and  then  run  PERCEPT.  Any 
time  the  " 

LOCATE  16,  6 

PRINT  "simulator  waiting  for  a  response  to  a  prompt,  you  can 
print  out  a" 
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LOCATE  17,  6 

PRINT  "screen  dump  to  an  EPSON  FX,MX  compatible  printer. 
Depress  " 

LOCATE  18,  6 

PRINT  "<SHIFT-PrtSc>  for  medium  resolution  (110  dpi)  or 
<SHIFT-PrtSc-h>  for" 

LOCATE  19,  6 

PRINT  "high  resolution  (240  dpi).  After  the  printer 
finishes,  answer  the" 

LOCATE  20,  6 

PRINT  "prompt  and  the  program  will  continue." 

LOCATE  23,  6 

INPUT  "Press  <RETURN>  start  the  simulation...",  A$ 

CLS 
END  IF 
END  SUB 


SUB  readin  (q, DATAFILE,  vect(2),  class(l),  n,  trainf ile$ ,testf ile$) 


OPEN  trainf ile$  FOR  INPUT  AS  #1 
FOR  I  =  1  TO  n 

INPUT  #1,  vect (I ,  1),  vect (I ,  2) 
INPUT  #1,  class (I) 

NEXT  I 
CLOSE  #1 

END  SUB 


SUB  readin2(q, DATAFILE, vect (2) ,class(l) ,o , trainf ile$ ,testfile$) 
OPEN  testfileS  FOR  INPUT  AS  #3 
FOR  I  =  1  TO  o 

INPUT  #3,  vect (1,1),  vect (I, 2) 

INPUT  #3,  class (I) 

NEXT  I 
CL0SE#3 


END  SUB 
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Experiment  1/2  Flatplate  17.88cm 


snf"  -  Table[(flbslsaa[li]]]''2  +  Rbs[sballi]]]~2  +RbsIsabIIUll/'2  +  Rbslsbbltil]] 


ListPlot[span,PlotJoined  ->  True,  P!otLabel->"SpanM] 
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Experiment  1/2  Dihedral  8.944cm  Flatplate  17.88cm 


span  =  Table[(Rbs[saa[[i]]]A2  ♦  flbs[sbaimiJ-2  +Rbs[sab[[i]]r2 


ListPlotlspan.PlotJoined  ->  True,  PIotLabei->"Span"] 


+  Hbs[sbb[[i]J] 
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Experiment  1/2  Dihedral  8.944  cm 


■  «r 


span  -  Table[(flbs[saa[[i]]]A2  +  flb$lsballi]]]*2  +flbslsablli]]]~2  +  flbslsbblli]]] 


ListPlot[span,PlotJoined  ->  True,  PlotLabel->*Span*J 
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Experiment  1/3  Dihedral  8.944  cm 


span  »  Table[(flbs[saa[[iJIJA2  ♦  flbs[sba[[i]]]-2  +Rb$[sab[[i]]]-2 


ListPlotfspan.PIotJoined  ->  True,  PlotLabel->'Span"l 


+  Rbs[sbb[[i]]l 
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Experiment  1/3  Flatplate  17.88cm 


span  -  Table[(fibs[saa[nil]A2  +  Abs(sba[[i]]]A2  +Bbs[sablli]]]>'2  +  flbslsbblli]]] 


ListPlot[span,PlotJoined  ->  True,  PIotLabel->’Span"] 


-Graphics- 
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Experiment  3d  Dihedral  89.44cm  Flatplate  37.888cm 


span  -  Table[(flbs[saallill]~2  +  flbs[sba[lilllA2  +Hbs[$abl[i]]]~2  +  flbsfsbblli]]] 


ListPlot[span,PlotJoined  ->  True,  PlotLabel->"Span“] 
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Experiment  3g  Dihedral  4.472cm  Flatplate  37.88cm 


span  -  Table[(flbs[saa[[i]]]~2  +  flbslsballi]]]A2  +flbslsab[[i]]]''2  +  Rbslsbbllill] 


ListPlot[span,PlotJoined  ->  True,  PfotLebel->“Span“I 
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Experiment  27  Dihedral  110cm 


span  =  TabIe[(Rbslsaal[i]]]A2  ♦  flbslsballi]]]"2  +Rbs[sab[Ii]]]~2  +  Rbs[sb b[[i]]] 


UstlPlot[span,PlotJoined  ->  True,  PIotLabel->*Span’l 
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Experiment  41  Dihedral  180cm 


span  =  Table[(flbslsaa[[il]lA2  +  flbs[sba[[i]]]A2  +flbs[sabl[ill]~2  *  flbs(sbb[[i]]] 


ListPlot[span,PIotJoined  ->  True,  PlotLabel-^Span"! 
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Experiment  86  Dihedral  80cm  Flatplate  20cm 


copole2Long=Table [2f [Ra  [q2c [ [i] 1 ] — 0, 

If  [Iza[q2c  [  [i]  ]  ]—0, 

N [ArcTan [-1 . 0] , 5] , 

-Sign  [Iatq2c [ [i] ] 1 ] *N [Pi/2, 5] ] , 

N [ArcTan [Re [q2c[[i]]],-Ha[q2c[[i]]]]  ],5  ] 

,  {i, pulses}] ; 

The  following  section  plots  the  co  and  cross  polaiizadons  on  the  Poincare  Sphere. 

copoleList=Table [ 

{copolelLat [ [i] ] ,  copolelLong[ [i] ] } 

, {i, pulses}] ; 

The  following  is  the  sphere  graphic 


test=Graphics3D [sphere] ; 
Show  [test.  Boxed.  ->  False]; 
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Experiment  91  Dihedral  100cm  Flatplate  20cm 


span  -  Tab!e[(Rbs[saa[[i]]]~2  +  flbslsbal[il]]~2  +Rbs[sabllil]]"2  +  Rbslsbb[[i]]] 


ListPlot[span,PlotJomed  ->  True,  PlotLabel->"Span"] 
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Experiment  95  Dihedral  120cm  Flatplate  10cm 


span  »  Table[(flbslsaa[[i]]]~2  +  flbsIsba[[i]]lA2  +Rbs[sab[[i]]]A2  +  Rbslsbblli]]] 


ListPlot[span,PlotJoined  ->  True,  PIotLabel->"Span"l 
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Experiment  99  Dihedral  120cm  Flatplate  80cm 


span  »  Table[(Rbs[saa[[i]]]A2  +  Abs[$ba[[i]]]A2  +flbslsab[li]]]~2  +  Bbs[sbb[[i]]] 


ListPlotlspan.PlotJoined  ->  True,  PlotLabel->"Span"] 
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Experiment  103  Dihedral  160cm  Flatplate  GOcm 


span  -  Table[(flbslsaa[[il]]~2  +  nbs[sba[Ul]l~2  +flbs[sab[Ml]-2 


Lis  tPfotfspan, Plot  Joined  ->  True,  PlotLabel->"Span"J 


+  fibslsbblli]]] 
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